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ABSTRACT 


A linearized mathematical model.is developed which simu- 
lates the dynamic response of the Navy F/A-18 for the study 
of flight control reconfiguration. Тһе аттстасс ΝΥ... 
as a multi-input multi-output, sampled data, closed system, 
which couples the dynamics of the flight control Systemes 
the aircraft linearized small perturbation egua bh The 
discrete time, state variable equations for the system are 
then formulated. A computer program is developed which will 
compose the model matricies and compute the response of the 
aircraít to stick and rudder input. 

To study flight control reconfiguration, the “model ape 
individual actuation of either a left or right control sur- 
face. Aircraft response to the actuation loss of either the 
left or right stabilator is simulated in the program p 
program is designed to implement the reconfigurable control 
mixer, Cürrentiy under e for the Self-Repairing Digital 
Flight Control System. 

The computer simulation was written in VS οκ ИНИ 
copy of the program and simulation results are included in 


the appendices. 
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I^. INTRODUGUTON 


Tactical aircraft face airborne and ground based threats 
which continue to grow in number and capability.” Іп Тир 
combat engagements the lethality of these systems will pose 
a considerable threat to ο. survivability. In аас 
to combat losses, the loss of aircraft in battle damage re- 
pair, or awaiting repair, will significantly redüce our tac- 
tical forces. Projections on the Survivability оик 
forces [Ref. 1] during the initial days of engagement 1mkEES 
cate 682 of the tactical aircraft will be Ogio action 
after the third day of battle. Of this total, 22% wm 
lost in combat and 462 lost while in repair or awam ы 
pair. This is indicated in Fig. 1.1 which comes from Kefer= 
ence 1. It is clear from projections such аз 1815, 1 а 
Porces ame to remain a superior threat to the enemy; οσοι E 
ued emphasis must be given to reduce the combat vulnerability, 
and increase the reliability and maintainability of οὐ. 
Пса ТЕСНЕЕ» 

To improve the combat effectiveness of future tace moe 
aircraft, the United States Air Force initiated R Om 
Repairing Flight Control System, Reliability and Мазза 
ability Program. Reference 2 outlines the program plan and 
goals. An Air Force sponsored study [Ref. 1] showed that 


Significant improvements in aircraft survivability and 
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Figure1.1 Tectical Aircraft Battle Damage Repair Statistics 


pot 


reliability can be achieved sn self repairing digital 
ТИРЕ control systemi 

One of the principle features of a self repairing system 
is flight control reconfiguration. A reconfigurable кеи 
control system automatically counters aircraft loss of Gone 
trol due to impairment, or loss, of a control effector mM 
fundamental idea is to utilize the existing undamaged effec- 
tors to preserve the normal flying qualities of the unimpaired 
aircraft. The approach outlined in Reference l is to design 
a reconfigurable control mixer to be placed between the flight 
control laws and the control effectors. In the face of air- 
craft impairment the mixer would reallocate control commande 
to the unimpaired effectors so that flight critical раш шы 
yaw, and rolling moments would be preserved to the maximum 
extent possible. Using the control mixer concept, the exist: 
ing flight control laws would πο ЕВЕ ЕЕ. 

Although the development of the self repairing system is 
intended for the advanced tactical aircraft, it is possible 
that the mixer could be implemented in existing airframes. 
e could be done with a control reconfíiguratuion modd 
interposed between the aircraft flight control computer cua 
fly-by-wire actuators [Ref. 1]. 

For this thesis the McDonnellF/A-I18 was chosen tovs 
the reconfigurable control law concept utilizing the contes ИИ 
mixer as described above. This choice was made based on the 


following assumptions given in Reference 1 for implementation 
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ОШОО ШО ο υπ τα τοῃ солсгої law process in an existing 
airframe: 
1) “The effector complement provides redundant effector 
systems and surplus control power for each flight 


ООШ ЛЇК оле! force and moment." 


oe the flient control is a full authority, fly-by-wire 
deco Ελ πρι control system." 


eye A control law design exists which has been carried 
Оев опе twnimparred airplane, but which is suf- 
πο οποιο febust that only first order impairment- 
induced stability derivative changes need be 
aecounted Ғос im a drop-in reconfigurable mixer." 
Based upon these assumptions it is felt that the F/A-18 is 
ΙΙ. Suited for the reconfigurable control law study. 

In this thesis a linearized mathematical model was de- 
veloped which simulates the dynamic response of the F/A-18 
aircraft to stick and rudder inputs. The model includes the 
ШО laws for the longitudinal, lateral, and directional 
axes for the cruise phase of flight. The fly-by-wire actua- 
tors and sensor n s r modeled and joined to the air- 
frame linearized, small perturbation model. The perturbation 
foe’ was obtained from the Simulation Control and Technology 
Ш сир, Flight Systems Branch, Strike Aircraft Test Direc- 
ШКО ас the Naval Air Test Center. A computer program 
was developed which, with existing software at the Naval 
Mt eraduate School, will simulate F/A-18 longitudinal, lat- 
eel, and directional response to stick and rudder TOUS. 

Тһе model was designed to implement the control mixer 


meen Concept outlined above. Future work at NPS will utilize 


the simulation program to develop algorithms for determining 


ΙΙ. 


статс stabiliey sand geome ccs 





Πιν ОЕ METHODOLOGY 


τον 

This chapter discusses the methods and assumptions used 
NEN rmustatevcawemathematdcal model which simulates the dynam- 
EE пзе of the F/A-18 aircraft. The flight control sys- 
Memes described, including the simplifying assumptions used 
to develop the control system model. An overview of the com- 
plete system, which couples the control system model with 
mie airtrame small perturbation model, is then given with a 


ШЕШЕСІ description of each functional component. 


ШЕ FLIGHT CONTROL SYSTEM DESCRIPTION 

ЕЕ deccripetomeonr the F/A-18 control system and 
Mmm@eory Of operation can be found in the Flight Control Sys- 
tem Design Report by McDonnell Aircraft Company [Ref. 3]. 
The following discussion briefly describes the basics of the 
mereht control system and the control law mechanization, and 
is intended to facilitate understanding the model development. 

ite primary flight control system in the F/A-18 is a 
fly-by-wire, full authority, control augmentation system. 
ο σοπετοι law computations are performed by four flight 
EN! computers operating in parallel. Each computer re- 
eet es input from the aircraft motion sensors, air data 
memputer, and pilot stick commands. The computer operates 


EN "c Input signals according to the control law algorithms 


η» 


апа μις the command signals to fly-by-wire electrohy- 
draulic servoactuators. Figure 2.1, taken from Reference 3, 
shows a functional block diagram of the flight control o e 
Exclusive of angle of attack and air data sensors, the system 
has quadraplex redundancy. The system provides two fail op- 
erate performance for augmented motion feedback control. A 
third failure causes the system to revert to either open loop 
direct electrical link control, or Stabilator πε ыш 

Соп го, 

The control augmentation system is gain scheduled with 
angele of attack and a@ir data tO provide optimum flying quali- 
ties throughout the flight envelope. Cross axis intercon- 
nects (e.g., rolling surface to rudder intereonneee) eam 
provided for turn coordination and maneuverability at high 
angles of attack. The control system also provides feedback 
to Counter nmerpbrat coupling at high roll rare к 

Тһе F/A-18 has ten primary flight control surfaces: 

Right and left stabilators, leading edge [lop p > MU 
edge flaps, ailerons, and rudders. Longitudinal control аа 
provided by collective stabilator, and collectives kc ESS 
and trailing edge flaps.  Lateral-directional control is 
provided by differential stabilator, differential сао 
and trailing edge flaps, ailerons, and NEUE E Collective 
leading and trailing edge flap deflections are scheduled by 
the control laws and are a function of angle of attache 


flap positions are designed to provide optimum DR Ka 
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Outer Loop Reterance Inputs trom MUX Bus 
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Functional Block Diagram of Flight Control System 


cruise, and improve flying characteristics in maneuvering 
and high angle of attack flight. Figure 2.2) (copied sie 
Reference 3, shows the control surface positions and direc- 
tion for positive deflect iong 

The objective of the thesis did not require a complete 
model of the F/A-18 control system as given in Reference 3. 
The following assumptions were made to reduce the model 
complexity. 

1) The aircraft is operating in the ир апа анау τα 
phase. Under this assumption the control laws are 
operating in the auto-flaps-up contiguration dM. 
trol law configurations for the takeoff and landing 
phases were not modeled. 

2) Only inner loop control is modeled. For inner) ioe 
control the pilot provides commands to the system. 
The auto functions (outer loop control v" r Kim 
modeled. 

3) Control is provided by the control augmentation Systeme 
The unaugmented modes such as direct electrical link 
or mechanical backup are not modeled. 

4) The failure logic provided to reconfigure the control 
laws in the event of a sensor or actuator failure is 


not modeled: 


5) The aircraft is operating with gear up, speedbrakes 
in, and no external stores. 


6) Spin mode control logic. is nor nod ο 

7) The aircraft trim system is not που μα. 

8) High angle of attack conditions are not considered in 
the thesis model, therefore control law configurations 
for this flight condition are not modeled. τ... 


thesis high angle of attack is defined as flight con- 
ditions above 15 degrees. 
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С. F/A-18 DYNAMIC ΜΟΡΕΡ ΟΥ, ΠΡ; 

The ЕЕ block diagram of the F/A-18 model which 
couples the flight control system to the basic airframe is 
shown in Figure 2.3. Thediagram represents a multi-input 
multi-output, sampled data, closed loop control system. 
Theory on the analysis of sampled data systems, as the one 
shown i Fig. 2.3, is extensive and covered in a пишБе ш 
texts (see references). The development of the F/A-18 model 
assumes the reader has a rudimentary understanding of control 
uheory and ın particular the theory of sampled data control 
systems. 

In the nomenclature used to represent the control signals 
pilot inputs are prefixed by the letter 'P'. Ассчиасвіпр 5 
nals, and signals produced by the aircraft sensors are pre- 
fixed by the letter 'E'. Control surface deflections same 
prefixed by the letter 'D'. The nomenclature used to repre- 
sent the aircraft perturbed motion variables and солее 
surfaces is given in Table 2.1. To denote the motion axis 
which is being controlled, the signal will be suffixed EE 
either x, y, or z to denote longitudinal, lateral, Or m ы 
tional axis respectfully. Finally a matrix or a vector яя! 
be denoted by an upper case letter. A scalar will be dénoted 


by a lower case letter. 


The input vector shown inm rio n E 


P(t)* = ірх(Е) py B SH v 
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БЕП V/A το Παρετα προτα (Luoton εἰς ӘТПБІД 


2] 


(3) =ХН= (3) 53 





(A)A (3) 53 (х) 53 


(3) 49+ (3) 5Х4= (1) $Х 
SOIWWNAG ЧО$МЯ5 а/ч 











(аҙаағ(а)шхн-(а)А Η 


(3)ез| 0 (хҙез (3) 23 
2 


(3) *XH-(3)qd (х) па+ (х) 2х9= (у) 23 
(3ЈА (3)α OU d (3)d 
(3) Q5 (3) ux as (3) ux (3) *39* (3) eXa- (3) ex 


OU n8 (X) >ХЧ= (Т+х) ЭХ 
SOIWVNAG S3WWH.hHIV SOIWWNAG HOLV(LOV 


v/a ΝΙ 59 «ΜΟΙ 'IONLNOO а/ч 


ТАРЫН! 


ΜΟΤΙΟΝ УАКТАВІЕ5В АМ) СОКТКОГ ЭПККАСЕ οἱ... 


4 Ditech wae ВЕБ, ЗЕ 
nz normal acceleration G 

aa angle of attack DEG 
yr yaw rate DEG/SEC 
ГГ roll rae DER SEC 
ny lateral acceleration G 
str right staba Paden. DEG 
stl left stabilator Q 
ter right trailing ed@err Ban $ 
tel left trailing edge flap n 


j right leading edge flap 
lel left leading edge flap 


аг righc ail ron Ú 
а1 ett ea ome. i 
ES right rudder š 
Т! lett rudder τ 


represents the aircraft longitudinal and lateral Stick jae 


rudder deflection in inches. Ihe output. cerene 
E 
Y(t) = [q(t) nz(t) aa(t) yr(t) rra ελ] а 


represents the perturbed motion of the aircraft about some 
Steady state operating condition. The motion variable mm 
are degrees, degrees/sec, and G's. Each block in the dia- 
gram contains a mathematical model which simulates the dynam-— 
ics of that particular component. The control 3 bM я 
contains the aircraft flight control law algorithms mode а 
as linear, time invariant, discrete state теи аа One 
processing channel of the flight control computer пезсті 
above is represented in the control law model. The input 
vectors to the control law model are the discrete stick ama 


rudder input. signals in Tine CS 
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t N 
DES БУО рск] 2030 
and the discrete motion feedback signals from the sensors, 


πμ πο ENERO yr(Ck)UEr(k) ny(k) ] 2T) 


(The units of the motion feedback signals are degrees, degrees/ 
БЕС апай G's). The input vector to the control law equations 


is therefore 
E 
вк) | Ce 
Mime Output vector from the control law block, 


век)" = [estr(k) estl(k) eler(k) elel(k) eter(k) 


etel(k) ear(k) eal(k) err(k) ὁ) 


Шергесепіс іПһе discrete command signals to the flight con- 
Πο! actuators in degrees. The actuator TT signals en- 
ter the GAIN block which represents the configuration gain 
matrix discussed in the introduction. The state space equa- 


mrmons Shown in the actuator block model the dynamics of the 


EE Cnr control actuators. The input vector to the actuator 
block, 
18 
DIM coc estbhot)veler(t) elel(t) eter(t) 
И саты! ealt) етті(г) егіс)! up 


mepresents the continuous time, actuator command signals in 


σσ The output vector from the actuator block 


p 


represents the control surface deflections in degrees, 


D(t)" = [dstr(t) dstl(t) dler(t) dlel(t) dter(t) 


dtel(t) dar(t) dal(t) drrot mM MI (278 


The deflection vector is input to the airframe small pertur- 
bation model represented by the state variable equations in 
the airframe block. The output vector from the small per- 
turbation model, Y(t), enters the sensor dynamics block ИИ 
contains the state variable model for the aircraft rate gyros, 
accelerometers, and angle of attack sensors. As discussed 
above the sensors output the feedback signals which are sent 
to the control laws. In Chapter III the mathematical models 
which simulate the dynamics of each component in Fic. M 

will be developed in detail. 


Analog to digital converters are modeled as impulse sam- 


plers. It is assumed that all samplers are operating at the 
same, constant sampling rate. (The actual system uses multi- 
rate sampling. In the thesis model only a single sample 
rate is used. The program actually allows any desired sam- 


pling rate to be input.) The mathematical operation per- 
formed by the impulse samplers is shown in Figure 2.4. 
Modeling the analog to digital converters as impulse samplers 
is a valid assumption if the quantization error of the ace 
system is at an acceptable level. Digital to analog conver- 
ters are modeled as zero order hold devices. The mathemati- 


cal operation performed by the zero order hold is shown in 
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ШИ О ОО ου ο ЕІпаііу ІС і5 аввишей that the processing time be- 
tween the sampler inputs and zero order hold output is very 
much less than the sampling period (i.e., no processing delay 
pime is assumed). This assumption is used in transforming 
the continuous time state equations (actuators, airframe, and 


sensors) into discrete time equations. 
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Figure 2.4 Mathematical Operation Performed by Impulse Sampler 
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Figure 2.5 Mathematical Operation Performed by Zero Order Heme 
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ЕТТЕ ТЕСЕБРОРМЕНІ 


СОНО ИЕЛИ МОВЕ DEVELOPMENT 


ΠΠ ΞΞΘΤΠΙ, 200 drrectlonal control law 


models were developed from the block diagrams and information 


Provided in Reference 3. In addition to the simplifying as- 


 ШПріісопс listed in Chapter Ii, many of the components in 


this system were eliminated based on the following 


memsiderations: 


1) 


2) 


3) 


4) 


5) 


6) 


The control law model will be coupled to the small 
perturbation model of the F/A-18. It is assumed 
that operation of the control system will remain 
within a linear region. Therefore the non-linear 
components in the control laws were eliminated. 
These include position limiters, rate limiters, 
and dead band regions. Other non-linear functions 
in the system which were essential to the model 
(pitch stick gradient for example) were linearized 
by a taylor series expansion. 


For the same reasons discussed above, the portions 
of the system which provide inertial decoupling 
were not modeled. 


In the model, the control signals are input as dis- 
crete signals. Therefore stick and rudder dynamics 
are not modeled. 


Noise, which may be detrimental to control system 
eaor monce due to aliasing, is not introduced into 
the model. Therefore the anti-aliasing prefilters 
were not included in the model. 


me ОООО ОО Л ЕЕ deS еге пос included in the Е/А-18 
πο Tame model. ittercihore tLhewstructural notch fil- 
БЕС οἱ που included in the model. 


To prevent discontinuities in the signals, the con- 
PeolmlawonWetlize taders in portions of the system. 
κ ο πε παν ιτας Gould Occur during start up, fail- 
οτι παπι είοηπς. Since none of these condi- 
tions are included in the thesis model, faders have 
not been modeled. | 


7] 


7) The F/A-18 control system uses multi-rate sampling 
in the input and feedback paths (200 900 3mm Ға 
sampling rates are used). To develop a state var- 
iable model only a single sampling rate was con- 
sidered. Therefore the iteration averagers, used 
to mathematically combine two discrete signals of 
different sampling periods were not modeled. The 
simulation program allows any desired sampling rate 
to be input. For this thesis 80 hz was used as the 
sampling rate: 
Figures 3.1 and 3.2 show the simplified block diagrams 
of the longitudinal and lateral-directional cont: O 3 
The inputs are the discrete stick and rudder commands, and 
motion feedback signals. The output signals are the commands 
to be sent to the flight control actuators via the recon е 
urable gain matrix. Together Figs. 3.1 and 3.2 make up the 
control тато оек Еа то kunt oe и 
The blocks in the control law model represent two basic 
transfer functions: Function gains, and digital ΠΟ. 7 
Table 3.1 lists the notation used to represent the transfer 


functions. Function gains and digital filters are десети 1 


below. 


ЖАБ ЕД ый 


CONTROL LAW TRANSFER FUNCTION ОЛ ОА 


Prefix Transfer fun on 
ru Function sa m 
LEE Longitudinal digital ETE M 
τ, ο Lateral digital τι s 
ү Directional digital? meee. 
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ПЕ ОЕ ЕКО 

ШООГО ОЕ ПШ Ἱῃπ (пе control law diagrams perform the 
Шы еш gain scheduling described in Sec. II.B. The functions 
operate on the air data, angle of attack, and normal acceler- 
ation to compute the system gains. Table 3.2 lists the nota- 
tion used to represent the function inputs. The mathematical 
equations which derine ane ΠΕΤΕΡ o given an Refer- 
ence 3. Copies of the functions used in the control law 
model are given in я A. In the simulation program the 
metion gains are computed using steady state conditions for 


EU input values. 


TAPCEE T? 


FUNCTION INPUT VARIABLES 


ps lcu Ἔπεσε pressure BS Гг? 
qc Dynamic pressure les t 
ri Pressure ratio (ps /qc) ND 

nz Normal acceleration G 
alpha Angle of attack DEG 


2e Piece ens 


Lead-lag filters are used in the system to shape the 
Output response and provide ТТТ gain and phase margins. 
ШІ іІпСесгасгог із used in the forward loop of the longitudinal- 
System to provide zero steady state error between command and 
Meedback. The control system design report [ÍRef. 3] gives 
mee tilter’s continuous time transfer function. For the con- 
Exc law model, the digital filter coefficients were computed 


Meme the Tustin transform. All filters were modeled as first 


9} 


order systems with two numerator coefficients and a single 


denominator coefficient. For example 
H(Z) = РЭМ ЕЕЕ (3 "> 
LZ - РЭО 
is the model for filter P9. Appendix B gives the dio Ы 


filters seda ene model, and the method used to compute the 
coefficients. 
37 ‘State Space Mode s 
Standard control system analysis techniques меге ШИ 
to derive the state space models for longitudinal and lateral- 
directional control laws shown in Figs. 3.199209 ашы 
following sections outline the procedures muc mE 
qe boned dies p us s m 
The block diagram in Fig. 53.1 Γαπεε ο. 
input/output signal ilo pario: 
1) Pitch rate to collective sta bm knon 
2) Normal acceleration. to collective Smr ERE 
3) Longitudinal stick to collective stab mie, 
4) Angle of attack to collective leading edge flap 
5. Angie ofat a kro sonik να. trailing edge flap 
To obtain the individual path transfer functio nm p m sign a E 
are mathematically combined to give the following three Z- 


transform equations: 


Estx(Z) = H1(Z)Eq(Z) — H2(Z)Enz(Z) = H3(Z) ix (2 


Elex(Z) = Πα (σσ κ. (3... 


οι - Hs 2) haa Z ) (3.4) 


NN pPresssons for each transfer function are given in Ap- 
Memdix ©. To obtain a state space expression for longitudi- 
KW πετοὶ laws, the individual Z-transfer functions are 
first expressed in state variable form. The state variable 
equations are then combined according to equations 3.2 - 3.4. 
mots procedure is as follows: 


ве сте чее function for the pitch rate to collective 
Stabilator path can be expressed as: 


И о но ве) 


(Z-P9D) (Z-P2D) 
The numerator coefficients are functions of the system gains 
and filter coefficients. The roots in the denominator are 
the poles from Jill j ji Q dd РЕ) <Appendix C gives the de- 
Mailed expressions for the numerator coefficients). 
2) The state space representation of Eq. 3.5 is obtained 


using the parallel programming method outlined in 
Reference 4. 


xl(k+1) PO DERO ІЗІ. | 

-| |] | + ББ 9230 ) 
х2(К+2) ТЕЛЕ СЕС 1 
Xu CO 


ИО = ев | 


| 4553 ed K) EI. 
x2(k) 


Applying the same procedures to the remaining Z-transfer 


D ons in Bq: 3.2 results in similar expressions: 
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x3(k+1) POD ο κ. 1 

| |- | Mese (3.8) 
ха(к+2)0 = 0 ο ὃν. ОЛИ 1 

ХЕ) 


estx2(k yn mmn on өзе | nzst3 enz(k) (322 


x lee) 
for the normal acceleration path. H2( 4) e200 


x5(k+1) = POD x5(k) + 1 px(k) (3.10) 


estx9(k) z pxstl x5(k) = ОЕ (3 ЖИ 


for the longitudinal stick рак Ы 
Appendix C details the procedures used to com- 


pute the coefficients in the output equations. The nomena 


ture used to represent the coefficients in the output equation 


combines the notation of the input Signal апа υπερ 
control surface, followed by a number indicating the coeffi- 


cient's numerical order in the equation. Бог си 
GStl qs b σου. 


are the first, second, and third coefficients in the ορ)... 
rate to stabilator output equation (EG. 3.70) ы 
system of nomenclature the respective signal path of the co- 
efficient is easily identified. 

3) The state variable equations for HI(Z), H2(Z) απα WE 


are now combined accordine to) Eaves 


estx(k) = estxl(k) — estx2(k) esM (3.1 B 
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пе Есе. 21, 3.7 ana 3.11 gives 


XEK) 
x2(k) 
estx(k) - [ast1 qoum sti St -pxstl]Dx3(k) + 
x4(k) 
ӘСЕ) 
eq(k) 
вене] | же рх(К) Е 
| еп=2 (Кк) 


EIN ου ος ο Ἱαξείοπ. Equations 3.6, 3.6, and 3.10 can 


be combined to give 


х1(к+1) P9D lle) СЕ О 

ο P2D 250 EUST e 

E k-1)/[- ΡΟΡ τπτ +|[о [рх(к) 
х4(к+1) рэр xa I) ШІ? (с 

х5(К+1) PODH x5(k) Ono ο τ.) 


The state space representation of equation 3.2 is given by 
Mes. 3.13 and 3.14. These equations output the collective 
Stabilator command given the longitudinal stick and motion 
feedback inputs. Note that the motion feedback and longitu- 
Enal inputs have been separated. This facilitates coupling 
Ише control law equations to the aircraft equations to be 
developed later. 
ΠΝ. suilasr state variable equations are derived for the 
Mugen or eattackwto collective flap path transfer 
πη Коп Па)” ала П5(2): 
Moe =P LID ХОСЕ ғ 1 еаа(К) (3.15) 


s ОТ ШК ате τος) r'aale2 еаа(К) 5) 
Bor the AOA to collective leading edge path H4(Z): and 


9 


| 


x7(k+1) = P12D ХТ(Е στ... (3.10 


etex(k) » aatel x7(k) ср зара” а (3.1 P 


for the AOA to collective trailing еаре ρει, 


5) The state equations for collective stabilator, sama 
collective leading and trailing edge flaps are now 
combined to give the following state variable model 
for the longitudinal Comer oO ВЕ 


х1(к+1) POD x ke) 
x2(k+1) ГО χα, 
х3(К+1) P9D x3(k) 
xd bs P5D | x4(k)] + 
29120 € 319 POD x ee 
x6(k+1) ΙΡ Lio (Kk) 
x7(k+1)] р122)х7(ю)) 





LAU O С 
100 О 
ое ею 0 
ОСО ете κ) p< (3.19) 
0 0 Oleaa(k) 1 
ΦΟΝ 0 
ша? О | 
x LCi 
х2 
estx(k) Gstl ~qst2 пле пле л D: E. 0 x3(k) 
en. О О О О О aalel 0 Іше - 
etex(k) 0 0 0 0 0] О aatellix5 ο. 
хб С 
x 71859 Wy 
qst3-nzst3 0 eq(k) -pxst2 
| О О aale? ssl. 0 pone» (3,205 
О О ааге2Иеаа(К) О 


Equations 3.19 and 3.20 are the state space representarnom 
of the longitudinal control laws shown in the block атас η 
in Fig. 3.1. The equations represent the longitudine P 


trol law model which computes the collective stabilator 
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КООШ апа Ес сесте flap commands, given the longitudinal 
Stick and motion feedback inputs. 
b.  Lateral-Directional System 

The procedures outlined above are applied to the 
BENE 05-2 to obtain the state space model for the lateral- 
directional system. (To complete the discussion of the con- 
trol law model these procedures will be briefly described.) 
ll individual path transfer functions are first MIS by 
ВЕ пешастса11у combining the signals in Fig. 3.2 to give the 


moltlowing equations: 


Esty(Z) = -H6(Z)Err(Z) + H7(Z)Py(Z) + H8(Z)Pz(Z) (3.21) 
Eley(Z) = -H9(Z)Err(Z) + H10(Z)Py(z) 22) 
В) = HIIOZ)EPT(Z) 4 H32CZ)Py(z) (3.23) 

Ea(z) » -H13(Z)Err(Z) « H14(Z)Py(2) ж Н15(2)Р2(2) (3.24) 


D E CZ ZH OZ Err(Z)--* 


ЕЛІ ОШООО Ур 7 11720277002(27) (3.25) 


These transfer functions represent the following input/output 
signal.paths, numbered respectfully: 

КОШ oll rate to differential stabilator 

Е в tosustferential stabilator 

pee Rudder pedal to differential stabilator 

0) Roll rate to differential leading edge flap 

10) Lateral stick to differential leading edge flap 

Mie Roll rate to differential trailing edge flap 


Meee Lateral stick to differential: trailing edge flap 


Э? 


13) .Roll rate to ai c ron 

14) Lateral stick to ат теш 

15) Rudder pedal to aileron 

16) Тат гаге го таев 

17) ЕСІ! таге νο τους» 

18) Шағегаі1 ассе!Тегастош б (4 ег 
19). Lateral stick to ruddar 


20) Rudder pedal -towr ttdde: 


The expression for each transfer function and corresponding 
State equation are given in Appendix C. Note that the dif- 
ferential stabilator, ailerons, and rudder contain the ο. 
fer functions for the cross axis interconnects (e.g. GR 
and H15(Z) represent the rudder to rolling surface inte r ы 
песе, апа H19(Z) represents the rolling surface tomma 
interconnect.) 

The state equations for the lateral and direction- 
al control laws are given as: 
esty(k) О -rrst 
И O ry e 


etey(k) O -rrte 
ea(k) O -rra 


pyst pzst 
еуг(К) +ру1е 0 ру(К) (3.2 
cu pyte 0 pz(k) 
eny(k) pya  pza 


CLOSS 


for the lateral system ο πα 


х8(к+1) Y3D 71:169 B 
x9(k+1) Ya) x9(k) 

ХІОСЕТІ) Y | το... 
xd Gk = Y5D UD 
x12(k-1) Y5DIX 12 0 F 
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E MC δα O0 
C [yr O- Ny κο 
ШЕГЕ” DNO sce) 
ο O Olleny(k)} +11 Oj}pzck) 
OGOTO DE 
iow 
moO le) 
er(k) =|-угг1 per EE II pr р2г4 х10(К) |+ 
x | Ok 
СЕ) 
τ; 2469 
Е > ву еге (е) |; L р2г2] 15222) 
eny(k) paky > 
for the directional system. Note that the lateral system is 
of order zero. No filters were included in the model. 
ο. ο πρ πι ας Еве Control Law Models 


ГЕШ oea О алаи 3.20, and 3.26 - 3.28 are com- 


bined to give the 3-axis control law model. 


eq(k) 
enz(k) 
de] Wee 12 12x1 12x6 eaa(k) INDORE) 
HI DEN MN buck δε рук) (3.29) 

err(k) р2(К) 
епу(К) 

estx(k) eq(k) 

elex(k) enz(k) 

etex(k) О а ο ο елар) р (К, 

К) = Сс Kc(k) + Dfc eyr(k) + Вс py(k) ο ο 

eley(k) егг(к) p2 (k) 

есеу(к) епу( К) 

ea(k) 

er(k) 


l ions 3.29 and 3.30 are written in terms of the matrix 


E Ucients in Appendix D. These equations are represented 


Sur 


by the discrete state equations in the controls ee 
παν ο 
3. Variable басп i; 

As a final step in the control law model development, 
the variable gain matrix is introduced. Тһе command signals 
which are output from the control laws (Eqns. 3.29 апа - 
are distributed to the right and left actuators accordi mE 
the diagram in Fig. 3.3. For the шипіпраігей атшгсгатт η 
individual gains in Fig. 3.3 are set to unity. πε τοῦ E 


ing matrix gain equation represents the diagram in Fig. 3.3. 


'GAIN' 

esc [ О О = | 0 0 0 0 

estl(k) 1 О О 1 0 О О Ollestx(k) 
eler(k) 0 | 0 0 il 0 0 Ollelex(k) 
elel(k) О 1 О 0 -1 0 0 о егех ие) 
er or 0 0 il 0 0 -1 0 Ollestv(k) 
etel(k) О О [| О О | О Olleley(k) 
еаг(К) О О О О О 0 -1 Oletey(k) 
eal(k) 0 0 0 0 0 0 i O flea(k) 
егг( К) О О О О О О 0 lMer(k) 
erl(k) 0 0 0 0 0 0 0 1 


(3. 315 
The matrix gain equation vill be recomputed for the impaired 


dere Gat πα 


Б. ADTRCRAFT MODELS DEVELGE b I | 

A 3-axis control law model has been developed which will 
operate on the stick and rudder inputs and the motion feed- 
back signals. The model outputs the control signals which 
are distributed to the aircraft actuators through (һе το... 


ble gain matrix. The next step in buildine the FyAm o τ᾿... 
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Figure 3.3 Control Law Command Signal Distribution and Gains 
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is to develop the state ЕКІ Ml equations for the ο. 
actuators, and sensors. 
1.  Airframe Model 

The airframe model, obtained from the Flight Systems 
Branch at NATC, Iis а linearized sm m пето раа model for 
both the longitudinal and lateral-directional modes. The 
dynamic stability and control derivatives were generated from 
MANTO SP AS ποπ. РасКаге. 


The aircraft state [O DE 


model is trimmed, unaccelerated l-g flight. The airframe 


State equations are given as: 


LONGITUDINAL 
u(t) χε ΙΕ, 4x3 dsc а, 
WCE) = Fx w(t) + Gx dlex(t) (3.328) 
q(t) ГЕ) dtex(t) 
thed(t) thed(t) 
q(t) 3x4 u(t) 3x3 dstx(t) 
Mace) EI NEU + Ох dlex ma @" 
aa(t) ην, dtes CEN 
thed(t) 
LATERAL DERECI LOM 

dsty(t) 
СЕ) A A NN Ce) 4х5 dley(t) 
ΤΙ, =r με + Gyz duc NN (3.3 P 
Diu DIOE да С) 
phi CO) phie) de e) 

dsty Gey 
Е В. в) Зеба ате, 
p(t) = КИ ο πα + Dyz dtey(t) (3.28 
ny(t) p(t) da(t) 
phi D) Ч) 


Airframe variable definitions and units are lisc oe 


3.3. The definitions of the stability and control deriva τν 
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Ἱ ΕΕ ο 3 


RESPIHRAMECVARTABLE DEFINITIONS 


iotation Variable Hii ES 

u(t) ШО аала velocity perturbation τε. 
w(t) perturbed normal velocity IP E 
G L) perturbed pitch rate Гас 65 
Е) perturbed lateral velocity Pts 
Ir (t) aireralt perturbed yaw rate abus 
BE t) aircraft perturbed roll rate rad/s 
ihed(t) ον a eper turbed i pitch angle rad 
phi(t) da rait perturbed roll angle i 
fee x( t ) τε, саш от deflection D 
dlex(t) eollectuive beading edpe flap deflection ~ 
dtex(t) collective trailing edge flap deflection " 
ast y(t) teeth hal srapimekor der lection 
dley(t) ues uale csdu edoe deflection d 
dtey(t) differential trailing edge deflection i 
da(t) η! БЕТЕ оа š 

r (t) воет α5ξεῖσει τοῃ к 

пе Сс) monum l acceleration ο ου 
Басе) angle of attack rad 
ny(t) lateral acceleration b ES 


and associated units, which make up the matricies in Eqs. 
EN through 3.35 are given in Appendix E. Note that the 
[mets Of the perturbation model are not compatible with the 
 Пасгогс ог sensors. (The actuator output units are in 
degrees, and the sensor input units are in degrees/sec, de- 
Mees, and G's.) The input and output variables for the 
Eturbation model were scaled in the simulation program to 
Meoperly interface the models. Combining equations 3.32 - 


3.35 gives 
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L. u(t) dstx 
w(t) 4х4 4х4 w(t) AE “хо 411 55 
qu Fx 0 СЕ. Gx 0 di ex 
thed(t) = --------- thed(t) + --------- dsty (37318 
VCE) 4x4 424 αν 4x3 Axon `. 
Ее 0 νο πι, 0 Gyz  dtey 
p(t) p(t) da 
РИШ С) p EU dr 
μπε, dstx 
aCe 3x4 Эха Е. Ох 3x5@ dlex 
De Hx O GM Dx 0 dtex 
aalt) = --------- thed(t) + --------- dsty (3.998 
r(t) 3x4 3A πα, Ira 3x5 dley 
p(t) 0 πρ o 0 Dyz dtey 
ny E) ; a ee) da 
πι κ, απ 


To study reconfigurable flight controls the aire a 
model should be capable of using the full set of control sur- 
faces available to produce the required forces and moments. 
To achieve this the control surfaces are split into indepen- 
dent right and left hand complements (i.e., right elevato 
left elevator, etc.). The equations are then coupled so that 
a complement of control surfaces used either collectively, 
differentially, or as a single side, will produce the apps 
priate moments. For example the stabilators deflected col- 
lectively will produce a pitching moment, deflected 
differentially will produce primarily a rolling попепс, M 
Single side deflected will produce, to some degree, moments 
about all three axes. 

The F/A-18 airframe modeled in equations 3.36 and 


3.37 inherently offers control coupling through the всаз ы 


tor, leading edge flap, and trailing edge flap Sur m F 
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Е олат ПОИ Ора could be achieved with the 
EBerous and rudders. At πο 5 Е was written 
control derivatives were not available on the longitudinal 
Peeeects of the ailerons or rudders. (The rudder is capable 
ШЕ гое-іп ог meee our and will effect the longitudinal re- 
Euse of the aircraft. This feature is normally used during 
takeoff and landing.) 

р ЕЕ ООСО surtace" detlections into right 
and left hand complements, the following equations are used 
СП compute the deflection inputs to Eqs. 3.36 and 3.37. 
Eo refer to Fig. 2.2 which shows the control surface posi- 


tions and corresponding positive deflections.) 


LONGITUDINAL DEECECTITONS СЕКА DIRECTIONAL LT DEFLECTIONS 
wk =. (4561 + dstr )/2 Coa el E Sr 
dlex = (dlel + dler )/2 dley = -dlel + dler 
dtex = (dtel + dter )/2 dtey = dtel - dter 

Е СОА 

ПЕ УСТА dgar)/2 
Mere +t and I correspond to right and left surfaces. These 


ιν.  1ο0π5 are rewritten in the following matrix format: 
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| 


dstx 
dlex 
dtex 


dsty 
dley 
dtey 
da 
ΠΤ 


| 


ο οσο ον 


| 


ο το ασε 


се 


OOOOH 


! 


ο = 


t 


OLOT © 


LONCE [n axe mS 
0 0 0 0 0 0 

5 0 0 0 0 0 
(ES 5 0 0 0 
ΙΑΤΡΟ distribution maei 
0 0 0 0 0 0 

= О 0 0 0 0 

Om. pal ak 0 0 0 

0 0 0 = 1 0 

0 0 0 0 Oey 5 


и om 


ες» 


detr 
ss | 
dler 
dlel 
dter 
dtel 
dar 
dal 
drr 


ατα 


dstr 
dst1 
dler 
dlel 
dter 
dtel 
dar 
dal 
dr r 
du 


(328% 


(3.995 


Note these equations are for the unimpaired aircraft only! 


Damage to one or more of the control surfaces will change the 


LONG and LATD matricies directly. 


Replacing 


wo tix σε ο ος: 


of equations 3.5599 9md я 


the Input vectors ив 


ing modified airframe equations: 


u(t) 
w(t) 
q b) 


ft hed (tb > 


κ 
DE 
p(t) 


ha Ct) 


Fim 
4x4 | 4x4 
Га 0 
4х4 | 4х4 
Оо Г г 


и (Е) 

ντ, ax 
η. Сх 
thed(t) + 
να) 4х5 
yr(t) Gyz 
πμ 

μμ. 
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Е" 


πας. 


3.36 and 3 8 


39 gives the follow- 


dstr 
dst1 
dler 
4161 
dter 
dtel 
dar 
dal 
dE 
dci 


(3.405 


Hm Dm as tr 


ШСЕ) dstl 
E t) 3x4 | 3x4 w(t) ES Ld] ou 
mz« t) Hx | О πίε, Dx LONG dlel 
aa(t) = --------- thed(t) + --------- ШЕСІ 
mE t) 3x4 | 3x4 v(t) о O Fi el (3.41) 
NC t) 0 | τν t» бугт ID gar j 
ny(t) wet t dal 
pim E drr 
drl 


ВЕ names of the individual matricies appear above the 
equations. 
ОО Actuator Model 

вес оез опе Гот the flight control actuators 
were given in Reference 3 and are listed in Appendix F. 

These transfer functions are low order approximations of the 
more complicated actuator models presented in Reference 3. 
The models were used in the F/A-18 rigid body stability anal- 
ysis described in Reference 3, and approximate the frequency 
@m@eracteristics of the higher order models out to 5 hz. 

ШОО ОШ р се Ene actuators into the modified air- 
meme model, the transfer functions are first put into state 
ШОШО ЫЫ е form. This procedure is outlined in Appendix F. 

The following equations represent the actuator state variable 


mode]: 


24х1 24x24 24x1 24х10 10х1 


б ы --ο3 Ea(t.) 2) 
ТОТ ШЛА 2a] 
ИСЕ) = На τ) ση 


EN 


мһеге, 


Ea(t)* - [estr(t) estl(t) eler(t) elel(t) eter(t) etel(t) 


ear(t) eal(€t) еггі 2 Ἀπ“ 
апа, 


D(t)? = [dstr(t) dsti(t) dler(t) dlel(t) dter(t) dtel(t) 


dar t EE drr(t) druü Ελ] 


The input vector, Ea(t), represents the input signal in de- 
grees from the control lavs, via the GAIN matrix. The output 
vector, D(t), represents the right and left control sur iai 
deflections in degrees. Equations 3.40 - 3.43 are combined 


to give the following equations of the airframe plus actuator 





modei: 
Be Xp(t) Gp 
8х1 8х8 | 8х10 10х24 8х1 8х10 
Xm(t) Fm | Gm Ha Xm(t) 0 1054 
----- = ----------------- ----- 4 ----  Ea(t) 
Axl хваро Да 24х10 (з 
Кл E 0 | Ға Да Са 
Нр Xp(t) 
q(t) 
Е) 8х1 
1153 6x8 | 6x12 12x24  Xm(t) 
ЛЕ = Hm | ра Ha ----- (3.4 
Dep 2:5 sei 
ПЕ = Е) 
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F. before the names of the matricies appear above the equa- 
τους. Тһе атага рі ое actuator model inputs the command 
E nalis to we actuators, and outputs aircraft motion. Since 
the control surface deflections in the perturbation model are 
in radians, it was necessary to scale the Hm and Dm matricies 
to interface with the actuator model which outputs deflections 
in degrees. This is done in the simulation program. 
E Sensor Model 

In the final T. of the development of the aircraft 
model, the state variable model for the aircraft sensors is 
incorporated into the airframe plus actuator model. Тһе 
sensor transfer functions are given in Appendix G. In state 


variable form the sensor equations are given as: 


τα 
πε) 
ІШЕ! TEx EI ПО I osa E. 55η 
Xs(t) = Fs Xs(t) + Gs met) 
p(t) 
ny(t) 
ος) 
enz ον 
еаа( ғ) ον νι 
eyr(t) = Hs Xs E) (СЭ ) 
ӘРІГЕ? 
eny(t) 


Mee sensor model inputs the aircraft motion variables in 
degrees/second, degrees, or G's, and outputs the correspond- 
MmemstesnalsS to the control laws in the same units. Combin- 
ЕС 545. 3.44 —- 3.47 gives the folloving model of airframe 


plus actuators plus sensors: 
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Saal 32x32 | 32x1 e 32x10 
ρα) DN 0 Xp(t) Gp 10х1 
----- т ----------------- ----- + -----  Ea(t) (3024 
11х1 11х6 6х32| ПО ІШУ НЕ 
ο σς Нр | Ер Xs(t) 0 
Hps Xps(t) 
eq(t) Ах 
env t) Ou | Бери) 
eaa(t) = 0 | Hs ----- (3.49 
eyr Ce) (ІЛЕСЕ! 
ορ ο) | XS ED 
πμ. 


The aircraft model inputs the actuator signals гош Ерс ΙΕ 
trol laws via the GAIN matrix and outputs the motion signage 


from the sensors which are sent to the control law qua GD ú 


С. ASSEMBLING ТИЕ O VARA TLZSYO [BHD I D 

A mathematical model for each component in the contre" 
system shown in Fig. 2.3 has now been developed. Before the 
individual components of the model can be assembled, the dis- 
crete state equations for aircraft model (Eqs. 3.48 and 3.49) 
must be computed. Performing this operation thema r c m 


state equations for the aircraft are given as: 


zx ü ena I E ТО 
Хрз(К+1) = Арѕ Kps(k) + Врз  Ea(k) (3.508 
6x1 6x44 aut 
FS k) = “Новы (3.5 
Where 
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Xps(k)* 2» [Хш(К) | Ха(к) | Хз(К) ] 


вк) S festr(k) estl(k) eler(k) elel(k) eter(k) etel(k) 


ear(k) eal(k) егг(К) ег1(К)] 


Meck)" = [eq(M) enz(k) eaa(h) eyr(k) err(k) eny(k)] 


The Aps and Bps discrete matricies are computed as follows: 


Aps =. ш о) 
ο i 
Bps = £ ος “ds X Gps 59.) 


Where ts represents the svstem sampling time. 

The GAIN matrix is now introduced to interface the dis- 
Eure aircraft equations with the control law equations. Re- 
Placing the input vector, Ea(k), in the aircraft equation 


ЕП the r.h.s. of the GAIN equation (Eq. 3.31) gives the 


В ом по: 
43x1 а b xÍ ο ο po] 
Хрз(К+1) = Aps Xps(k) + pps ОВ ЕСК) о 
6x1 6x44 ВЕ 
K k) = Hps Xps(k) 3 5 


For convenience the control law equations are repeated 


below: 
ie x 1 Jo E ey eom Gl Бора 
КС(Е-І) = АСЕР хс К) ce Bfc Es(k) + ο. ερ, ΕΟΝ 
8х1 ЛЕЛ στι 8x6 6х1 τη πε] 
Шк) = Ce πας) τα κα οι κ) + МЕЗ ЗБар 30) 


ο] 


Мһеге: 


Хс(к) с = [Xx(k)|Xz(k)] 

Ез(к)° = [eq(h) enz(k) eaa(h) eyr(k) О (КОСО О 

P(k)" = [рх(к) ру(к) рг(к)] 

Ec(k)" » [estx(k) elex(k) etex(k) esty(k) eley(k) etey(k) 


еа) бегот m 


Equations 3.52, 3.53, 3.29, and 3.30 are now compe ВИ 


give the following matrix equation: 





43x1 43x43 43x10 10x8 8x6 6x43 ;, 43x10 10x38 е MI 
Xps(k+1) Aps + Bps GAIN Dfc Hps | Bps GAIN Cc Xps(k) 
12x1 12x6 6x43 | 12х12 1250 
ΕΠΙ Вне | АС Хс(К) 

a ox SNO ον 
pps GAIN ο . 
+  --------------- Б (3.50 
T233 pack) 
ВС 
6x1 6х43 | 6x12 43x1 6x3 
Еѕ (к) Hps | 0 Xps(k) Ор 
------ ш ----------------- ------ + --- py(k) (373258 
8х1 5:0 ο δν ον и 8x3 pz(k) 
БО К) Die Sips | Ce Xek) De 


These equations model the dynamic response of the F/A-18 


System Shown “in Figure 2... 


D. MODELING ЕРЕБСПОК ТИЕ |І 
Effector impairment is divided into four groups termed 
‘effector impairment classes’ (EIC) [κ6Ε. ΤῚ» Ῥ ΓΡ 3.19 


copied from Reference 1, defines the EIC Sand sin ee peer 
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CONTROL STABILITY & 
AFFECTED COBTRSCISAFPRECTED 











ACTUATION LOSS, 
SURFACE FLOATING 


ACTUATION LOSS, 
SUREAGS GEMIEPED 


ει ο 


PARTIAL SURFACE 
Γ᾽ 


PIqure 3.4 Effector Impairment Classes 


a 


their affect on aircraft stability and Control Shee eee 
Fig. 3.4, effector impairment class one is unique in t 
only aircraft control is affected. mihien e impair neH 
classes affect aircraft stability as well as "eon. comme? 
self repairing control system must be capable of detecting 
and classifying effector damage to compute the proper recon- 
ЕО Е. 

For the case of EIC=l, one or more of the aircraft con- 
trol derivatives will be altered. In the F/A-18 model de- 
veloped above this is reflected in the LONG and LATD matraeaee 
(Eqs. 3.38 and 3.39). For example if the right stabila NEN 
impaired the elements LONG(1,1) and LATD(1,1) would δε οσο... 


Zero: 


LONG 
0 2 0 0 0 0 0 0 0 0 
0 "ges 5 0 0 0 0 0 0 
0 0 0 0 ο а 0 0 0 0 
LATD 
0 | О О 0 0 0 0 0 0 
0 0 1 -1 0 0 0 0 0 0 
0 0 0 το) ll 0 0 0 0 
0 0 0 0 0 Dus | 1 О 0 
0 0 0 0 0 0 0 0 D" MES 


Compare the above matricies with Eqs. 3.38 and 3.39 fore 
undamaged aircraft. Note that the damaged system will now 


produce a lateral input for a given longitudinal command and 
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mice versa. For the model developed ieee sere sass only 


elass one effector impairments are considered. 


DE CONCLUSION 

In this chapter the mathematical models for each compo- 
nent in the block diagram of Fig. 2.3 were developed. The 
Ma vidual models were then assembled to form the complete 
model of the F/A-18 dynamical system. Next the simulation 
program is developed to compute the model matricies in Eqs. 
3.54 and 3.55, and compute the response of the system to 


wc k and rudder inputs. 


553 


IV. PROGRAM DEVELOPMENT AND МОРЕ ολ... 


A. INTRODUCTION 

To validate the F/A-18 system model a computer program 
was written to compose the model matricies in Eqs. 3.54 and 
3.55, and compute the system response to stick and rudder 
inputs. In addition: the proc M simulates an actuation loss 
of the right or left stabilator. The program was written in 
VS Fortran on the IBM 3033 computer at the Naval Postgraduate 
School. The program is organized to offer flexibility Ын 


future development and modification. 


Be PROGRAM ο ο 5 
The program may be divided into four major operations: 

Ly Data input 

2) πο απο ο с ошро ано 

3) Composition of the system matricies 

4) System response computation 
The operations are performed by a series of subroutines as 
shown in the flow diagram in Fig. 4.1. Tn σπα ЖБ 
solid lines indicate control flow, and the dotted lines um 
cate data flow. The subroutines perform the steps outlined 
in Chapter II to compose the system matricies and compute 
the response. A functional description of each sur n EE 
shown in Fig. 4.1 and associated variables is given in Appen- 
"τν (Шеше subroutines at NPS were used to perform the 
required matrix manipulations [Ref. 5]. These subroutines 


are also defined in Appendix H. 


E 
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Figure 4.1 Simulation Program Flow Diagram 


oy 


Execution of the program is controlled by the EIS NP 
program. The Exec ‘eee the input and output data files as 
shown in Table 4.1. Non eee options are provided to define 
data input/output files. The result of each subroutine соп- 
putation is written to the FA18 RESULT file СИ 
viewed by the user after program execution. A copy of the 


Exec program is given in Appendix I. 


TAREE а 


INPUT/OUTPUT FILE DEFINITIONS 


Record Device File Type Use 
Number Name 
01 Disk FAIS DATA Contains the ҒА18 program 


input data 


Do Disk ΡΑΙ 8 ΠΕ ΠΕΙ Contains the resulti 
the operations performed 
by each subroutine 


03 Disk OPTMATD DATA Contains FAIS8 system 
matricies (Eqs. 3.549 6 
3.55) formated οι... а 
control system design 
package at NPS. 


O4 Disk OP Teer DATA Contains the system 
response data computed 
by the FA18 program 
formated for the inter- 


active plotting promi m 
at NPS. 


Ἱ Паца пр 
Subroutines 'FLITE1' and 'FLITE2' perform m eee 


input operations. All data is read from the FAIS8 data file 
except the control and failure parameters which are ор 


interactively. 
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-- 


ΙΡ; en te aircraft steady state flight 
@emditions and control and failure parameters. As indicated 
ЖЕНІС. 2.1 the flight conditions are read from the ΕΔ18 
ВАТА file. The control and failure parameters are read in 
Eusractively. This allows the user to conveniently run the 
program for various control inputs and control surface fail- 
Mees. "'FLITE2' reads in the basic airframe stability and 
Entrol derivatives from the FA18 DATA file. The stability 
and control derivatives are arranged in the matrix format 
Іп in Eqs. 3.33 - 3.36, and in Appendix E. 

2. Air Data Computations 

Prior to composing the system matricies, air data 
computations are performed by subroutine 'AIRDAT'. The sub- 
ШОО Пе computes the air data inputs to the control law func- 
tions using the standard atmosphere equations [Ref. 6]. The 
Eu boram does not include the logic to compute atmospheric 
conditions above the gradient (troposphere) region.  There- 
Ere the computations are valid only up to 36,000 feet. 

3. Composing the System Matricies 

As seen in Fig. 4.1 the system matricies are com- 
posed by ten subroutine operations. The operations and 
associated equations in Chapter III are as follows: 

опросе the LONG and LATD matricies Eqs. 3.38 & 3.39 


ШІ Сопрове (һе modified airframe matricies Eqs. 3.40 & 
2 


ΙΙ. onpose Lhe actuator matricies Eqs. 3.42 & 3.43 


2 


4) Compose the sensor matricies Eds: 9 25 ҒИ 


5) Compose the airframe plus actuator matricies Eqs. 
Μπ; 


6) Compose the airframe plus actuator sensor matricies 
Eqs. 3.48 & 3.49 


7) Compute the discrete system matricies Eqs. 3.5008 a 

8) Compose the GAIN πα τι ος n N 

9) Compose thecontrol law matricies ας. ου 

10) Compose the total system maricies Eqs. 3.54 & 1 
The basic component matricies are composed by first generat- 
ing a null matrix of proper dimensions, and then assie Tomis 
the amc πο О PINO to the proper elements. These matri- 
cies include: 

1) ΟΝΕ ἀαπα ILADD 

D uc 

3) Sensor 

a Controle tay 

>) GAIN 
The coefficient values to actuator and sensor models are 
written into the program as constants; they are not contained 
on a separate data file. The coefficients in the LONG, LATD, 
control law, and GAIN matricies are first computed gnaw 
assigned to the appropriate matrix element: 

The LONG and LATD matricies are composed by Sub 

tine 'LONLAT'. As explained in Section III.C these mac momin 
reflect control surface damage for one of ІГпе ТОП: СӘРЕ 


surface impairment classes. In the present program only 
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class one effector impairments are simulated. For EIC=1 the 
mppropriate control strface coefficients in the 'LONG' and 
'LATD' matricies are set to zero. ΠΠ ОШ О ат п< the logic 
ENUnpair the right of left stabilator.  '"LONLAT' also com- 
poses the unimpaired airframe control matrix, GmO. This ma- 
ο will be used to compute the impaired gain matrix in 
Embroutine 'VGAIN'. 
The control law matricies are composed by subroutine 

ИИ’. The subroutine may be divided into four operations: 

ШО Compute function gains 

E Compute the filter coefficients 

See Compute the matrix coefficients 


4) Assign coefficients to the control law matricies as 
πολ σος ο ὅ υ.5. 


ШПЕ ратпс аге соприсеай ассогдіпр to the function definitions 
given in Appendix A. Тһе functions use the air data computed 
EN TRDAT', steady state AOA, and normal acceleration to 
compute the gains. All functions are programmed exactly as 
shown in Appendix A except for the following non-linear 


EUDnctions: 


IMU ee 1 ποπ ο Pome tewdinale Stick, gradient 
ов в ео. Lateral stick gradient 
Pine P Tos la Rudder pedal gradient 
КОШО onmesg RSRI non-linear gradient 


These functions were programmed nem dRe linear terms of a 
ШИЕ От Series Expansion about the origin. Some of the func- 
tions computed in the program are not used in the control law 
model. They are provided for user information on system per- 
formance. These functions are: 
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Function 24 Trailing edge flap schedule 
Punmerroma ° Trailing edge flap schedule qc limit 


Runot ron® у Leading edge flap schedule 

Punet ions 29 Leading edge flap schedule ri limit 
Pune cron o Nz limit on AOA feedback 

F tre EST Ї Rolling surface limit schedule 


Function 112 Lateral acceleration aan 
Function 113 Lateral acceleration gain 


The filter coefficients were computed using the procedures 
outlined in Appendix B to transform the analog Γπτιτστο M 
digital filters. The control law coefficients were computed 
using the procedures πο. in Appendix C for transformesis 
the control path transfer functions into state variable ГИ 
The function and coefficient values are output to the FA18 
ΗΠ file. 

Subroutine 'VGAIN' composes the GAIN matrix based on 
control effector impairment. For the unimpaired system ο 
GAIN matrix appears as in Eq. 3.31. The subroutine is 86-2 
signed to implement the reconfiguration algorithm described 
in Reference l. 

Prior to composing the total system matrici1es ame 
discrete form of the continuous System matricies Ποπ... 
computed. This is done by subroutine "EXPINT Jl ROSE 


"EXPINT' computes the matrix exponential, 


Fps*ts 
з р 


and the integral, 


eFPS"Sgg 


Qu cr 
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The discrete system equations arethen computed as: 


c QPPS*ts 
and 
ts % 
Bps = ге Р "ds x Gps 


fae Aps and Bps matricies are written to the OPTMATD DATA 
memes The data file is formatted for the control system 


design package at NPS. 


А. Response Computations 


System response is computed for the recursive equa- 
tion 


xs αἱ. οσο 55x1 ο ο ο 
ΠΠ 5 κ ο. + BF18 U(k) 


Ns response is computed for 500 data points. The AF18 and 
ο matricies are defined in Eq. 3.54. The state and input 


vectors are defined as: 
X(h)* » [ Xn(k) | Xa(k) | XsCk) | Xc(k) ] 
E 
Εμ 


Response data for all 55 states is written into the 
OPTPLOT DATA file which is formatted for the interactive 
EE Ing routines at NP5. The variables contained in the 


state vector which are relevant to the user for viewing are: 
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= Ер x2d =з 
ΛΞ хо = рт x20 cM 
м: х= сик x25 = dar 
хо UMS κ ο n x27 = ай 
NE = в x29 = агг 
xo a Хх10-- (161 χ31 eM 


The output equation (Eq. 3.55) was πο ποτ me 


С. PROGRAM TESTING AND MODEL VALIDATION 
To test the simulation program, four SetS ОК ΗΤΟΙ 
made for-the following cont oB par x 
1) Positive longitudinal stick, ao fami 
20) Positive longitudinal stick, rieshc srabilar т ы ын 
3) Positive lateral stick? uo гае 
4) Positive rudder pedal,-no Газ № es 
All deflections were .l inch step inputs of 3 second du Aee 
All runs were made at .6 mach/10000 feet. The sampling rate 
was set to 80 hz, (Additional runs were made а hz sam- 
pling rates with no noticable difference in the output re- 
sponse.) Aircraft steady state data, including the stability 
and control derivatives, are given in Appendix Ju p ll p om 
RESULTS file, including the function and concro Eme MEN 
efficient values, is given in Appen PEN T 
The model was verified for: 


1) Correct direction of motion of Cont mo lsum as аа 
corresponding aircraft mot nonk 


2) Expected aircraft response for a rues 
failure, 


3) Ргорег augmented aircraft om s m 
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Reference 7 contains response plots for the aircraft for simi- 
MN UIsHt conditions. These plots were used to verify the 
meet for the proper augmented motion. Some of the plots 
were reproduced for the thesis and are given in Appendix L. 
Mme control inputs used in Reference / are slightly different 
than those used for thesis model. Some comparisons however 
ean still NS All response plots for the thesis model 
Eu given in Appendix M. 
j ТЛ г response со а positive longitudinal input 

ШЕГІ по failure condition is shown in Figs. M-1 - M.3. 
ime motion of the right stabilator, and leading and trailing 
edge flaps compare favorably to the McDonnell model response 
и in Pie. Lei. Іп the thesis model the right stabilator 
initially travels -2.5 degrees. After a .5 second transient 
period the stabilator continues to -3.6 degrees in 3 seconds. 
Mme corresponding McDonnell response shows an initial travel 
of -2.2 degrees. After the .5 second transient period the 
stabilator continues to -1.8 degrees in 3 seconds. The 
major difference between the two responses is the stabilator 
rate, the thesis model being slightly greater than the 
Mewonnell model. (Possible explanations for the discrepancies 
in the model are given below.) 

The leading and trailing edge flaps, driven by angle of 
ο τσκ feedback, deflect approximately 7.5 degrees in 3.4 
Ends as shown in Figs. M.1 and M.2. The corresponding 


McDonnell response show the leading and trailing edge flaps 
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to deflect approximately ^ degrees in 3.5 seconds. (Note the 
initial conditions of the flaps show Fig. PP lares s ΟΡΟ 
as in Figs. M.1 and M.2.) Aircraft pitch rate anum MM 
angle response for the thesis model are shown in Figs. M.2 
and M.3. The thesis model achieved a maximum pitch rate of 
8 degrees/second compared to 4 degrees/second for the 

MOTOREN STI model. Also the McDonnell model achieved a con- 
stant pitch rate in 2.5 seconds. The thesis model achieved 
a 2.5 degree pitch attitude in 5 seconds compared to 15 
degrees shown in.Fig. 1.1. The increased trave MES EER 
stabilator shown in Fig. M.1 (compare Fig. L.l) may cx p P WE 
the discrepancies in the pitch and flap responses. 

Figure M.4 shows the response of the aircraft with the 
right stabilator failed. Note the decrease in pitch rate 
and pitch angle magnitudes. Figure M.4 shows a negative 
roll rate which is to be expected for a right τοπ ο. 
failure. Reference 7 did not give aircraft response for 
GQOntrol surtace failures. 

Aircraft response to a positive lateral stick Τη. 
shown in Figs. M.5 - M.8. The motion of the ο απο. 
aileron, and trailing edge flap compare favorably to the 
McDonnell response shown in Fig. L.2. The adileron n ШЕ 
thesis model achieves a steady state deflection of 1.1 
degrees in .5 seconds. The McDonnell response shows a s eci 
state deflection of 1.6 degrees in .5 seconds. Тһе (БЕ 


model stabilator (Fig. M.5) deflects 0.25 ἄθρποςς Πα. 
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seconds, the DUIS LEE McDonnell response shows a 0.6 
lu codem lectionain 0-5 seconds. Тһе trailing edge flap 
EUN. M.6 deflects 0.4 degrees in 0.3 seconds which com- 
pares favorably to the response shown in Fig. L.2. (Leading 
edge differential flaps are not used at the Πα ο οοπαἳ Έτος 
ЖЕСІ Punction sich sets the gain on the leading edge 
MU path, is computed as zero.) The rudder response δια. 
M.6) shows mini ià negative deflection which is expected 
Bor coordinating a right turn. The magnitude of the deflec- 
tions for the two models is approximately the same, however 
the thesis response is slightly more oscillatory than the 
McDonnell response. The aircraft roll rate, yaw rate, and 
Ш апеле response are shown in Figs. M./ and M.8. Again 
only slight discrepancies exist between the two models. 

line aircraft response to a positive rudder pedal input 
[emonown in Figs. M.9 and M.10. The corresponding McDonnell 
response is shown in Fig. L.3. Comparison of the two models 
shows the shape of the responses to be approximately the 
same. However the magnitudes in the thesis model are very 
much less than the McDonnell model. For example the rudder 
achieves a maximum deflection of 0.07 degrees as shown in 
ШИЕ. 9. Compare this to a 1.5 degree deflection for the 
E unell model shown in Fig. L.3. The response of the 
differential aileron and stabilator for the thesis model 
Mees. M.9 and M.10) is essentially zero. This is expected 


mere the rudder to rolling surface interconnect gain, 


6 


func tror o т аа as zero for low angle on attacks 
(See Appendix A, function 39.) It is Unknown ee а 
and stabilators respond as shown in Fig. L.2 Ε555 5), 
response shown in Fig. L.3 is based on a different flight 
control program than that used in the thesis model. 

Possible explanations for the discrepancies απ MM a 
sponse between the two aircraft models are: 

1) Improper derivation of the simplified control law 
model based on the assumptions and procedures given 


in Chapters anchi. 


2) Programming discrepancies resulting in erroneous 
computations. 


3) Errors in the aircraft linear small рег и шаш ш 
model. 


4) Differences between the modeling techniques sae 
in Reference 7 and in this thesis. (1Ππ6 ΜΟΌΛΙΝ ΠΤ. 
included all aerodynamic and control system non- 
linearities, as well as the effects due to digital 
time delay and quantization.) 

5) The McDonnell model response made available to the 
author may be based on a slightly different tlio 
program than used in the thesis model. 

The author investigated each of the items listed above. It 
is felt the assumptions used in developing the control law 
model are correct given the available information on the 
F/A-18 control system. The simplifying assumptions used aim 
developing the model may account for some of the differences. 
The program code was thoroughly reviewed and revealed no dis- 
crepancies. It is assumed however that this item remains a 


possible source of error. Two possible sources of error 


exist with respect to the aircraft реггог ато s s 
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№) The coefficients in the NATC small perturbation 
model are computed with respect to the aircraft 
body axes (see Appendix E). This may effect the 
magnitude of the motion variables which are fed 
απο το theNeontrol! laws. 


2) The aircraft model developed in Reference 7 con- 
siders the offset position of the accelerometers 
ШЕСІ Petair ralt center of gravity. It has- not 
been determined if this effect is considered in 
the NATC perturbation model. 


rt her investigation into these possibilities should be 


Genducted. 
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ТУ. CONCLUSIONS AND RECOMMEND 11. 


A mathematical model was developed which simulates the 
dynamic response of the F/A-18 aircraft.: mme non Е 
designed to implement the Ше - gain matrix fae 
the study of reconfigurable control systems. A program 
was written which composes the discrete time state variable 
matrices or the атгегатеЕ "ande omues eS response to 
stick and rudder inputs. The P a sn also simulates the 
actuation loss of the right and left stabilators. кесро 
plots of the thesis model were displayed and compared to 
the aircraft model developed in Reference 7. Possible 
sources of error were discussed, and it is recommended that 
further investigation into each of these areas be conducted. 

τε is also recommended that the NATC program which com- 
putes the stability and control derivatives be acquired Eui 
the Aeronautical Engineering Department at the Naval Post- 
graduate School, to be made an integral part of this program. 


This would expand the ability of the program to simulate the 
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APPENDIX B 


DIGITAL FILTER MODELS 


ШЕ PITCH RATE LEAD-LAG FILTER P2 


ШИЕ 1-Е22,65 + 1 
DENN Nu Le (БО) 
мса! 


Using the tustin transform let 


2 (7-1) 
.ν. -τ- 
ts(Z+1) 
which gives 
(P2N1) ( P2N2) 
2(.015)(1+F22) + ts КЫ o (015) (1+P22) 
------------------- ЕЕ | 
2(.015) + ts ο δι τες 


2(.015) + ts 
ΤΕ ο) 
( P2D) 


BEC NORMAL ACCELERATION LAG FILTER P5 


P5(S) = --------- (B.3) 


Using the Tustin transtorm: 


(PSN1) (PSN2) 
ts ts 
See А σα = 
08 + ts 08 + ts 
Р5(7) = --------------------- (В. 49 
OB - ts 
[MEE MS ] 
OR + ts 
(PSD) 
Ge INTEGRATOR 589 
{| 
P9(S) = (απ 59 
5 


Reference 3 approximated this filter using the backward 


difference method. Let 


7-1 
S = ----- 
tsZ 
which gives 
(РОМ1) (Ρ9Ν2) 
ғе? + 0.0 
pou шы €—e——— —— m = — =~ = μ.ο 
7-41 

(pop) 
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ЕШ) LEADING EDGE AOA LAG FILTER Pll 


1 
НЕ... (B.8) 
δρ οος τι 
Шезпа the Tustin transform: 
(Р11М№1) (Ρ11Ν2) 
Es mS 
E ss TZ + тт 
"ΙΖ τ επ .39(2) + ts 
Р11(7) - ------------------------------ (В.9) 


.39(2) + ts 


(РІІГ) 


РШТЕАТГІКС ЕГСЕ FLAP AOA LAG FILTFR ΡΙ2 


Р12($) = --------- RATO) 


Using the tustin transform: 


(P12N1) (P12N2) 
ts ts 
ΕΞ ο ως 
.79(2) + +5 79(2) + ts 
Р12(2) - ---------------------------- (B.11) 
70(2) - ts 
στο ыы S= 


(P12N3) 
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Е. YAW RATE CANCED5F τ᾽ 


(Y3N] ) ( Y3N2) 
2 -2 
Болма с ln πε ες 
2 + ts 2 + ts 
Y3 (Z) = 22-2 
2 - Ев 
meo xc. Dm 
2 + ts 
(Y3D) 
6. RSRI LEAD LAG FILTER ΥΠ 
755 41 
ү5(5) - --------- 
5S + 1 


Using the Tustin transform: 
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25275) F tS ts =- 20 75) 

== == === ЕЕЕ 127 4 ----------- 

2С pr 2(.5) + ts 

ҮҰ5(7) ----------- 222. аа 
2(.5) - ts 
7 s S ee 5 
ο ο -- 

(YSN3) 


114 


(В. 128 


(В. 139 


(В.14) 


(В. 058 


APPENDIX C 


SIGNAL PATH TRANSFER FUNCTIONS AND STATE SPACE MODELS 


ШИ РІТСН RATE TO COLLECTIVE STABILATOR PATH H1(Z) 


From Fia. 3.1 the forward path transfer function is 


БСО О - [PPO*P?*F68*F32a*F12 4 Р2%((Е68%ҒЗ2а)ғҒ40)! (с.1) 


Let: 
А - F68*F32A*F12 
B = (F68*F32a)+F40 


Rewriting Eq. С.1 using the filter Z-transform expressions 


for P9 and P2 gives 


A*(PON1*7Z + PON2)*(P2N1*7 + P2N2) 


(Z-P9D)*(Z-P2D) 


B*(P2N1*Z + P2N2) 
ee Е Gr) 


Шарапсіпа Ес. С.2 


ЕШ (7) 


{fA*PON1*P2N1 + B*P2N1] Z 
+ fA*( PON] *P2N2+P9N2*P2N1]) + B*(P2N2-P2N1*POD)] Z 
+ fFA*PON2*P2N2 -— B*P2NO*POD]] / (Z-P9D)*(Z-P2D) 


беге? 


ΠΡ» 


Equation C.3 can be expressed in terms of the followinoc 


partial fraction expansion: 


2 
bO*Z + Біо 2 qstl acta? 
H1(Z) 9» ---------------- = gqst3 + ------- + ------ (С.4) 
(Z-P9D)*(Z-P2D) (Z-P9D) (Z-P2D) 
Where: 
ασε» - ου 
2 
qstl = (bO*POD + bl*(P9D) + b2) / (P9D-P2D) 
2 
qst2 = (bO*P2D + b1*(P2D) + b2) / (P2D-P9D) 


Equation C.4 can now be put in the followinq state variable 


form: 
х1(к+1) ΡΟ О іхі1(Ҡҡ) 1 
= ЦЕ (к) ( COEM 
x2(k+1) ο ο ο 1 
х1(К) 
estxl(k) = із ast2] + aqst3 q(k) ίσο 
x2 15) 


The sytem of nomenclature used for the coefficients in the 
output equation (Eq. C.6) is described η... τ᾽ 


In the remainder of the appendix onlv the first and last 
Steps of the state variable development will be given for 
each transfer function. The equations for each of the 
coefficients are given in the simulation proaram listina. 


то 


МӨБМАП АССЕТЕКАТТОМ ТО COLLECTIVE STABILATOR PATH H2(Z) 


БИІ - (рІ2%Ба82а%3.5%Р9%Р2 + 3.5*F32a*PS |] (ο и) 
х3(к-1) РЭР О x3(k) | 
= + nz(k) (ο 5) 
x4(k+1) O P5D x4(k) 1 
x3(k) 
estx2(k) = nzstl nzst2 Е ЕЗИК) КЕС 
x4(k) 


PITCH STICK TO COLLECTIVE STABILATOR PATH H3(Z) 
H3(Z) = [7.0*F32a*F12*P9 + F20*F32al ο 


Note:  F20 defines the longitudinal stick gradient. The 
value 7.0 is a taylor series approximation 
about the origin of F20. (See F20 in App. А) 


x5(k+1) = POD x5(k) 4 1 px(k) (Со) 


В — РБ ХК) + pxst2 pxik) Cord 
ANGLE OF ATTACK TO COLLECTIVE LEADING EDGE FLAP H4(Z) 
H4 (Z) = 1.328*P11 η 


Note: The leading edge flap schedule is defined by F27. 
ве тие еа Еа С. 13 is the l.e. flap 
gradient. (See F27 in App. А) 


x6 (k+1 ) 


PllD x6(k) + #1 aalk) (ССА) 
elex(k) = aalel x6(k) + aale2 aa(k) εδ eS) 


ΤΙ; 


ANGLE OF ATTACK TO COLLECTIVE TRAILING EDGE ЕГҒАР НИ 


Н5(2) = 12305 l (C. Tey 


Note: The trailing edge flap schedule is defined hy 
F24. The value used in Eq. C.16 is the t.e. flap 
gradient. (See F24 in App. A) 


x7(k+1) 


PI2D x7(k). * 1 аа(К) (С. 0 


etex(k) = aatel x?(k) ща а (COHEN 


ROLL RATE TO DIFFERENTIAL STABILATOR PATESHOZ 


H6(2) - rv7*(FA4-rk6t) (С. 1599 
estyl (k) = ву (Ако επ (C.20) 
where: 


rv7 - MIN(F6*F35,F6-F101) 
rk6t = LIMIT(0O,.12-F4,.12-F4) 
DE UL X 
LATERAL STICK TO DIFFERENTIAL STABILATOR PATH πι (72) 


H7CZ) m 3229* rv T Eb Ae ШИ (С. 21 


Note: The lateral stick aradient is defined bv Fl. 
The value 3.22 is a taylor series approximation 
about the origin of Fl. (See Fl in App. A) 


esty2(k) = 3.22*гу7*Е7* (ГЕЗ ТЕЛ НЙ (С. 2 


LES 


RUDDER PEDAL TO DIFFERENTIAL STABILATOR PATH H8(Z) 


H ΕΕ 33*rv7*F39 22297) 


Е ЕО КОЕ aed ose 7 V7 “FSS pz) (С.24) 


ROLL RATE TO DIFFERENTIAL LEAPING EDGE FLAP PATH H9(Z) 


HZ HOS ТЕЛ ЕСЕЛІ) 16725) 


eleyl(k) - F93*(F44rk6t) rr(k) (C26) 


LATERAL STICK TO DIFFERENTIAL LEADING EDGE FLAP H10(Z) 


H10(Z) = 3.22*F93*F7*(F13+F4) [5273 


eley2(k) = 3.22*F93*F7*(F31+F4) py(k) (c 28) 


ROLL RATE ТО DIFFERENTIAL TRAILING EDGE FLAP PATH H11(2) 


H11(Z) = F31*F34*(F4+rk6t) 25, 


etevl(k) = F31*F34*(F4+rk6t) rr(k) кесе, 


LATERAL STICK TO DIFFERENTIAL LEADING FDGE FLAP Н12(2) 


ҚАУ Е 3 22 ХЕ 7Ж<ЕЗІҰЕЗА4Х(ҒІ 3ҒҒЕА) ΓΕ 5 


etey2(k) = 3.22*F7*F31*F34*(F13+F4) py(k) (сез2) 
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ROLL RATE TO AILERON РАНО! 


μις ο Е35%Ғ36" 5 το, (C. 399 


eal(k) F35*F36* .5*(F44+rk6t) rr(k) (C. 332 


LATERAL STICK TO AILERON PATH H14(2) 


Н14(7) 3.22*P7*F35*5F36*-5* PME (С.З 


еа2(К) с κο ο 
RUDDER PEDAL TO AILERON PATH H15(Z) 

Н15(7) = F14*F39*1.33*F35*F36 (С.З 
еа3 (к) » F14*F39*1.33*F35*F36 pz(k) (c.3 

YAW RATE TO RUDDER PATH H16(2) 

H16(Z) = F45*F96*COS(alpha) *Y3 (C.39) 


Note: In the simulation program alpha is set to the 
steady state angle of attack in degrees. 


xB(k+)) = Y3D xP(k) r] vI P p 


erl(k) = vrrl x8(k) Усти 
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о. 


R. 


ROLL RATE TO RUILER FATH E17(2) 


H17(Z) = F45*F96*Y3*SIN(alpha) 


+ F45*F38*Y5*F30*0.76*((2*rra)+(2*rrst)) (С.40) 


Note: The second term on the r.h.s. of the eauation 
results from the rollino surface to rudder 
interconnect path. The value, 0.76, represents 
F42, the RSRI nonlinear gradient. It is a 
taylor series approximation about the origin 
of F42. 

x2 (k+l) pO О τος. 1 
= + pru (C.41) 
x3(k+1) 0 Y5D x3(k) l 
КаК) 
Smee) = rrr rrr + rrr3 rrlk) 7; 
x3(k) 


Ie 2) 


POTTERA ACCELERATION TO RUDDER PATH Н18(7) 


Edo E90 ο 5) 


er3(k) - nyr ny(k) C44) 


S. LATERAL STICK TO RUDDER FATT tae) 


H19(Z) = F45*Y5*F38*F30*0.76*((2*pya)+(2*pyst)) (С 


x4(k+1) 


er4(k) = ругі хак р а (GR 


Т. RUDDER PEDAL TO RUDDER PATH Н20(7) 


H20(Z) = F45*Y5*F38*0.76*F3O0*((2*pza E 


- F45*F14*( .5-(F17*F114) ) ία 


Note: The first term on the r.h.s. oft the egui iei 
results from the rudder to rolling surface 
interconnect path. 


x5 (k+1) YS5D хӨСЕ) к с Fp 


er5(k) = ρε. κο κο ο στο τα 


YSD X4 κ wu ΙΙ 2E (CR 


. 491) 


46) 


dM 


48). 





APPENDIX D 


CONTRƏL LAW MATRICIES 
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APPENDIX E 


ENADIPDITYCANDOCONTROL DERIVATIVE DEFINITIONS AND UNITS 


The matricies in the small perturbation model (Eas. 3.32 
through 3.36) are defined in this appendix along with the 
dimensions of the stability and control derivatives. This 
information was obtained from the flight systems branch at 


the Naval Air Test Center. 


FX MATRIX 
1 2 3 
ДІ XU XW XO-WR 
ZU ZW (ZO-UB) 
MM eee 
1-ZWD ] —- 7 WD 1-ZWD 
MWD* ZU MWD*ZW MWD* ( ZO-URB) 
3 MU-*------ MW+------ MQ+----------- 
EZ WI ] -ZWD IS 
4 0 0 Al 


A 


4 


ea CTET 


1-ZWD 


MWD*ZD 


S 


GX MATRIX 


HXEMATRIX 


ος 





DX MATRIX 


P YZ MATRIX 


ΠΗ; 


СЯН 


jo 


3 
О 
Z DRE 
l-2WD 
О 
5 4 
YP+WB G CHT 
l-YVD jc AVE 


NVD* (YP+WP ) NAD GG IPP 


GYZ MATRIX 


1 2 3 4 5 
YDHT YOLE XE YDA YDR 
1-YVD 1-YVD 1-YVD ] -YVD 1-ҮУГ 
NVD*YDHT NVD*YDLF NVD*YDTF NVD*YDA NVD*YDR 
NDHT4------- πη, πε, == ΝΡΑ!------- NpR+- aa 
1-YVD 1-YVD 1-YVD 1-YVD 1-YVD 
LVD*YDHT LVD*YDLF LVD*YDTF LVD*YDA LVD*YDR 
LDHT+------- LOLE- ва LDA+------ LDR+------ 
1-ҮҮГ 1-ҮҮГ 1-ҮУГ 1-ҮУГ 1-ҮҮГ 
0 0 0 0 0 


HYZ MATRIX 


| 2 3 4 
1 О | О О 
2 О О 1 О 

Yu YR YP 0 
: | со E 
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PYZ MATRIX 
Т 2 3 
О О О 
О О O 
YDHT YCE Y DTF 
А 1-YVD l-YVD 


d 5 

O O 

O O 
YDA ҮрЕ 
1-YVD 1-YVD 


The following abreviations were used in the above 


matricies: 


WB = 


CTET 


STHT 


BODY AXIS LONGTUDINAL WIND 


BODY AXIS VERTICAL WIND 


COS (THETA) 


СІМ(ТНЕТА) 


χυ 


XQ 


Х058 
XDTH 
XDS 
XDLF 
XDTF 


1058 
HDTH 
MDS 

МОЕ 
МОТЕ 


ШАУ) 
YVD 
YR 
YP 
YDA 
YDR 
YDLF 
YOHT 
YDTF 


DEFINITION OF D 5 ERIVATIVES 


R 


RIYATI 


LONGITUDINAL FORCE 
LONGITUDINAL FORCE 
LONGITUDINAL FORCE 


OF; 


LONGITUDINAL FORCE 


LONGITUDINAL FORCE 


LONGITUDINAL FORCE 
LONGITUDINAL FORCE 
LONGITUDINAL FORCE 
LONGITUDINAL FORCE 


VERTICAL FORCE 
VERTICAL FORCE 
VERTICAL FORCE 
VERTICAL FORCE 
VERTICAL FORCE 
VERTICAL FORCE 
VERTICAL FORCE 
VERTICAL FORCE 
VERTICAL FORCE 


PITCHING MOMENT 
PITCHING MOMENT 
PITCHING MOMENT 
PITCHING MOMENT 
PITCHING MOMENT 
PITCHING MOMENT 
PITCHING MOMENT 
PITCHING MOMENT 


PITCHING MOMENT 


LATERAL FORCE 
LATERAL FORCE 
LATERAL FORCE 
LATERAL FCRCE 
LATERAL FORCE 
LATERAL FORCE 
LATERAL FORCE 
LATERAL FORCE 
LATERAL FORCE 


ROLLING MOMENT 
ROLLING MOMENT 
ROLLING MOMENT 
ROLLING MOMENT 
ROLLING MOMENT 
ROLLING MOMENT 
ROLLING MOMENT 
ROLLING MOMENT 
ROLLING MOMENT 
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WITH RESPECT TO: 


FORWARD VELOCITY 
VERTICAL VELOCITY 


PITCH RATE 


VERTICAL ACCELERATION 
SPEED BRAKE 


THROTTLE 


HORIZONTAL STABILATOR 
LEADING EDGE FLAPS 
TRAILING EDG: FLAPS 


FORWARD VELOCITY 
VERTICAL VELOCITY 
PITCH RATE 

VERTICAL ACCELERATION 
SPEED BRAKE 


THROTTLE 


HORIZONTAL STABILIZER 
LEADING EDGE FLAPS 
TRAILING EDGE FLAPS 


FORWARD VELOCITY 
VERTICAL VELOCITY 


PITCH RATE 


VERTICAL ACCELERATION 
SPEED BRAKE 


THROTTLE 


HORIZONTAL STABILIZER 
LEADING EDGE FLAPS 
TRAILING EDGE FLAPS 


LATERAL VELOCITY 
LATERAL ACCELERATION 


YAW RATE 
ROLL RATE 
AILERON 
RUDDER 


DIFFERENTIAL LE FLAPS 
DIFFERENTIAL HORIZ STABS 
DIFFERENTIAL TE FLAPS 


LATERAL VELOCITY 
LATERAL ACCELERATION 


YAW RATE 
ROLL RATE 
AILERON 
RUDDER 


DIFFERENTIAL LE FLAPS 
DIFFERENTIAL HORIZ STABS ` 
DIFFERENTIAL TE FLAPS 


UNITS 


yoL* FT/SEC2 


Ε1/565᾽ 

RAD/SEC2 
FT/SE P 

RAD/SEC2 
РСТ/5Е 2 
RAD/SEC2 
RAD/SEC2 
RAD/SEC2 


ҒТ/ЗБ Ж 
FT/SE o 
RAD/SEC2 
ҒТ/ЗЕ А 
RAD/SEC2 
PCT/SEG2 
RAD/SEC2 
RAD/SEG2 
RAD/SEC2 


FI/SEGZ 
FT/SEGR@ 
RAD/SEC2 
FT/SECO 
RAD/SEC2 
PCT/SEGC 
RAD/SEC2 
RAD/SEC2 
RAD/SEC2 


ҒТ/5Е 2 

FI/SEC 

RAD/SEC2 
RAD/SEC2 
RAD/SEC2 
RAD/SEC2 
RAD/SEC2 
RAD/SEC2 
RAD/SEC2 


FT/SECZ 

FT/SEGS 

RAD/SEC2 
RAD/SEC2 
RAD/SEC2 
RAD/SEC2 
RAD/SEC2 
RAD/SEC2 
RAD/SEC2 


INITI OF DIMENSIONAL STABIL RIV Eo 


(Continued) 
МАМЕ PARTIAL DERIVATIVE OF: WITH RESPECT TO; UHITS 
NV YAWING MOMENT LATERAL VELOCITY FT/SEC2 
NVD YAWING MOMENT LATERAL ACCELERATION FT/SEC2 
NR YAWING MOMENT YAW RATE RAD/SEC2 
NP YAWING MOMENT ROLL RATE RAD/SEC2 
NDA YAWING MOMENT AILERON RAD/SEC2 
NDR YAWING MOMENT RUDDER RAD/SEC2 
NOLF YAWING MOMENT DIFFERENTIAL LE FLAPS RAD/SEC2 
NDHT YAWING MOMENT DIFFERENTIAL HORIZ STABS RAD/SEC2 
NDTF YAWING MOMENT DIFFERENTIAL TE FLAPS RAD/SEC2 


ЕИ! 


ДРРЕВЕТХ ЕВ 


ACTUATOR TRANSFER FUNCTIONS 


А. STABILATOR 


5 2 2(0.068) 


(----- ) + -------- ΕΤ 
Dst(S) 82.9 о 
Est(S) S ο ο τν) S 2 21ο 
[(----- ) + ------- S + 1G ) + ------- στα 
36.4 36.4 105.3 1.05 x2 


The procedures outlined in Ref. 8 were used to put the 
transfer functions in State variablen ornk 


1) Arrange the transfer function as 


2 
b2 S + b3 S + b4 
4 3 2 
S + al S + a2 S + a3 S + a4 


52 = ЭШЛЕ ОЗ ο πμ. 
Ыз προ * ao οσο. 
DA S IMMO ον α ο» 


al = 1.5410E+02 
2) Express the transfer function as a differential 


equation 
Dst(t) + al Dst({t) + a2 Dst(t) + a3 Dstlt) т ада по 


= b2 Est(t) + ЪЗ ЕБЕСЕ) Ба Ве 
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3) 


00 
ША) 


U3 
04 


x 
х2 
х3 
х4 


Ds 


Compute the state equation coefficients 


ко, 

= bl - al*UO = 
ey ee | 

= p3 — a]*U2 = 


b4 - al*U3 - 


О 


- a2* 


amo 
a2* 


UO 
Ul 
Ве 


2.377 E+O3 


a3*U0O = 


2520539250095 


иво = 2-7277E-07 


Arrange in the followinc state variable form 


Fs 
O 1 
= О О 
О О 
-a4 -a3 -a 
Bist 


ια I9 1059 


t 


xl 
х2 
х3 
х4 


αι. ο СО 


x 
х2 
x3 
x4 


Gst 


E 
+ 02 
03 
04 


БӘСІРЕ) 


Note input and output for all actuator models is in 


degrees. 


as above. 


The remaining state variable models are developeG 


5: 


В. “'LEADING EDGES SP I F 


Dle(S) J) 
Ele(S) S S 
(~--- + 1.0) ( ---- + 1.0) 
26.9 82.9 
Хо СЕ) I(E) 
= Fle + Gle Ele(t) 
МОЕ ко) 
ХИСЕ) 
Ге се) = Hle 
x21 E) 


С. TRAILING EDGE ELA? 


Dte(S) "ο 
Ete(S) S 2 ГОН: 
(----) % ------- 5 + 1.0 | 
35.0 3570 
| 
СЕЈ ХЕ) 
= Fte + Gte Ete(t) 
ХЕ” x2(t) 
ЕЕ) 
Dte(t) - Hte 
xou 


134 





Pe ALLERON 


Da(S) 1.0 
EalS) | S 2 2(0.59) 
се ре ------ о 
75.0 75.0 
RU E) x1(t) 
= Ба + Са Fa t) 
х2(%) КОСЕ 
СЕО 
Dalt) Ha 
xo B 


E. КОППЕК 


Dr(S) 10 
Er(S) 5 2(0.69) 
pcm --- S + 1.0 
yo poo sd 
πι 5 ІШСЕ) 
= Pr + Gr Erít) 
x2(t) πε) 
ΠΕ) 
Darte) = Hr 
XC) 
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The actuator state variable поде1 пағтісіев in no ΠΠ 


& 3.44 are arranged as follows: 


FA MATRIX 
Fst 
Fst 
Fle 0 
Fle 5 
Fte 
Fte 
0 Ба 
Ё Ға 
Fr 
ET 
GA MATRIX 
Gst 
Gst 
Gle 0 
Gle С 
Gte 
Gte 
Ga 
0 ба 
т Gr 
Gr 
НА МАТВТХ 
Hst 
Hst 
Hle 
Hte 0 
Hte M 
Hle 
Hle 
0 На 
A: Ha 
Hr 


НГ 
136 





APPENDIX G 


КОШЕ δὴ ο ӘК ТКАКӨКЕК ЕПУСІІОМӨ 


КОО RATE GYRO TRANSFER FUNCTION 


S 
----- + ITO 
Y(S) 31.27 
Е(5) $ $ 
[----- * 1.0] [----- + 1.0] 
167.8 461.7 


БЕ OCCELEROMETER TRANSFER FUNCTION 


S 
----- TL o0 
s (S) 235.8 
F(S) 5 2 2(.96) 
E CNN LLL μις 
305.3 Sob 


mee ANGLE OF ATTACK SENSOR TRANSFER FUNCTION 


Note: Transfer functions input dearees/second (or degrees) 
and output degrees/second (or degrees). The state variable 
model for the transfer functions in Eas. 3.46 & 3.47 are 
developed using the procedures outlined in appendix F. 


ре 


APPENDIX H 


SIMULATION PROGRAM SUEROULINGs 


Subroutines) “Pilar 


Reads in steady state flight conditions, 5... 
and rudder commands, and control summe 
failure parameters. 


Description: 


Calling sequence: CALL FLITE] (mach,alt,alpha,nz,ncont,nst, 


пер, апр. este Ен 


Input arguements: None 
Output arguements: 
mach Steady state mach number 
alt Steady state altitude in feet 
alpha Steady state AOA in dearees 
az Steady state normal acceleration. Normally 1.0. 
ts Sampling time in seconds 
neon Control nunhem: l=Longql tudinal stern 
2=Lateral stick 
3=Rudder deflection 
nst, nstp Control start time, and control вор ЕШ 
Nst and nstp are actually the iteration number: 
Start time (secs.)snst*ts 
Stop time (secs. )=nstp*ts 
amp Control amplitude given as deflection in inches 
roto CI Right and left stabilator failure parameters 
respectfully: No failure = О 
Failure = 1l 
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Eu5routine:  FLITE2 
Operation: Reads in basic airframe matricies from the 
ΤΙ: ile. 
EU Une arguments: CALL FLITE2(Fx,Gx,Hx,Dx,Fyz,Gyz,Hyz,Dyz, 
Nfx,Ngx,Nhx,Ndx,Nfyz,Ngyz,Nhyz,Ndyz) 


Input arguments: None 


eme put arguments: 


BE GX, Hx, Dx Longitudinal state matricies 
в, Сух, Нур, Dyz Lateral-directional state matricies 
Nfx,Ngx,Nhx,Ndx Two-dimensional vectors givina the 
ШЕУ?, Мау?, МҺһуг2, Мау? number of rows and columns of tne 
respective matricies. Fxample, 
Nfx(1) = Number rows in Fx 
Nfx(2) = Number columns in Fx 


199 


Subroutine: 


Description: 


LONLAT 


Generates the LONG and LATD matricies 
(Eqs. 3.38 and 3.39) based оп tne contre 
surface failure parameters. 


Computes the unimpaired input matrix, CmO, to 
the modified aiframe equations. The GmO matrix 
will be used in the reconfiguration algorithm. 


Calling Arguments: CALL LONLAT(rstf,lstf£,ifail,Gx,GyVZ ЧОЕ 


LATD,GmO,Nlong,Nlatd,NgmO,Nax,Nayz) 


Input arguments: 


тыс л ер 


ifail 


Gx GYZ 


Мах, στα 


Right and left stabilator failure parameters 
respectfully. (See subroutine FLITE1) 


Failure flag: No failure = 0 
Failure = 1 


Longitudinal and lateral-directional input 
matricies. 


Two dimensional vectors giving the number of 
rows and columns of the Gx, and, Gyz 
matricies. 


Output arguments: 


LONG, LATD 


GmO 


Nono, πιστα 
NgmO 


Matricies described in Section III.B. wb 
split the control surface deflections ο... 
right and left hand sides. The LONG and LATD 
matricies reflect control surface failure or 
damage. 


The unimpaired input matrix to the modified 


airframe equations (Fq. 3.40). Gm0 is 
composed using the unimpaired LONG and LATD 
matricies. It is suhsequently used in the 


reconfiguration algorithm. 


Two dimensional vectors givina the number of 
rows and columns in the respective matricies. | 
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Subroutine:  AIRDAT 


Description: Computes the air data inputs to the control 
laws based on standard atmosphere conditions. 


Calling arguments: CALL AIRDAT(mach,alt,temp,rho,psi,a, 


iab) 
Input arguements: 
mach Steady state mach number 
alt Steady state altitude in feet 
Output arguements: 
temp Standard atmosphere temperature in deqrees 
rankine 
rho Standard atmosphere densitv in 


ШБ=сес UM ft**A 


55 1 Standard atmosphere pressure in lb/ft**2 
a Sonic velocity in ft/sec 

ас Dynamic pressure in lb/ft**2 

ri Pressure ratio - qc/psi 
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Súübroutine: MOLEO 


Description: 


Composes the modified airframe matricies (Eos. 


3.40 апа ο 


Calling Arguments: 


CALL МОПЕО(Ғх,Сх,Нх,рх,Ғу?,Су?, НУ? ГЕ 


ГОМС, АТР, Ем, Сбт, Нм, Рм, МЕХ, МаАХ Т МАХ УСЕ 
Nfyz,Ngyz,Nhyz,Ndyz) 


Input arguments: 


Ex , Сх , ы / Ох 
πο ο. ο ουν 


LONG, LATD 


Nix Nox, NOX. Nex 
Nivz, Novz,. Nhyz, Wey z 
μπε»; Ni ads 


Output arguements: 


Fin; Gm, Hm, Dm 


Nim, Мат, Каш, Мат 


Longidudinal and lateral-directional 
basic airframe matricies (Бас и 
- 3.35) Read in from subprorr n аа 
FLIGHTI 


LONG and LATP matricies (fas. E 
and 3.39) Generated in subroutine 
LONLAT 


Two dimensional vectors giving the 
number of rows and columns іп ГЕ 


respective matrix. 


Modified basic airframe matricies 
(Eds. 3040 anos d 


Row and column vectors as described 
as described above 
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Subroutine: SENSOR 


Description:  Composes the sensor matricies (Eqs. 3.47 and 
3.46) 


Calling arguments: CALL SENSOR(Fs,Gs,Hs,Nfs,Ngs,Nhs) 
Input arguments: None 


Output arguements: 
Fs, Gs, Hs Sensor matricies 


Nfs, Ngs, Nhs Two dimensional row and column vectors 


Subroutine: ACTU 


Description:  Composes the actuator matricies (Eqs. 3.42 and 
C NES 


Calling arguements: CALL ACTU(Fa,Ga,Ha,Nfa,Nga,Nha) 
Enput arguments: None 


EUCDut arguments: 
Fa, Ga, Ha Actuator matricies 


Nfa, Nga, Nha Two dimensional Row and column vectors 
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Subroutine: PLANT 


Description: 


Composes the basic airframe plus actuator 


matricies (Eqs. 3.44 апа ЗЕ 


Calling arauments: 


Input arguments: 


Pm. Gm, Hm: ΠΠ 


Ба, Са, На 


МЕп, Ngm, Nhm, Ndm 
Nfa, Nga, Nha 
Output arguements: 


ое RD 
КЕр? Nep, Nie 


CALL PLANT(Fm,Gm,Hm,Dm, Fa,Ga,Ha,Fp,Gp,Hp, 
Nfm,Ngm, Nhm,Ndm,Nfa,Nga,Nha,Nfvo,Nap,Nhp) 


Modified basic airframe matres 
(Eas. 3.44 and 3.45) Genera я 
Subroutine: | ObEoe 


Actuator matricies (Еаѕ. 3.42 ала ЕВ 
Generated in subroutine ACTU. 


Two dimensional row and column σος. 


Airframe plus actuator matricies 


Row and column vectors 
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Subroutine: LAWS 


Beseription: Computes the function gains (App. A) 


Computes the coefficients in the control law 
matricies 


Composes the control law matricies 


Sling arguments: CALL LAWS(alpha,nz,psi,qc,ri,ts,Ac,Bfc, 
Bocce te, We. Nac, Note, Nec, Nec, Ndfc, Ndc) 


Duit arguments: 


alpha Steady state angle of attack in decrees 
n Steady state normal acceleration - 1.0 
psi standard atmosphere pressure at steadv state 


ее о Пе 


ас dynamic pressure in lb/ft**2 
n pressure ratio 
ts sampling time in seconds 


L ut arguments: 


Ac, Bfc, Bc απο) O тпасгяіслес (Есе, 3.29 апа 3.30) 
me, Dic, De 


Nac, Nbfc, Nbc Row and column vectors 
mee, Ndfc, Мас 
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Subroutine: 


pescerip on: 


VGAIN 


Computes the variable dain matrix based on 
control surface impairment. The subroutine 
15 disigned to use reconfiguration algorithm 
introduces Uu re E 


Calling arguments: CALL VGAIN(ifail,ifix,Gm,GmO,GAIN,Nom, 


NgmO, Ngain) 


Input arguments: 


ται] 


εινα 


GmO 


Ngm, NgmO 


Control surface failure flag set in subroutine 


-LONLAT ; No failure = О 


Falliro m 


Fix parameter. The GAIN matrix will be computed 
fo the impaired aircraft if τικ τα ane 

ifail = 1. Ifix is read interactively by sub- 
routine ELITEI. 


Airframe control matrix. From subroutine MORTER 
Unimpaired airframe control matrix. Generated 
in subroutine LONLAT. Note imao neo NE UID 


control surfaces are damaded or failed, Gm=Gm0. 


Row and column vectors 


Output Arguements: 


GAIN 


Ngain 


Variable gain matrix described in section 


ΕΕ ον 


Row and column vector 
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Bieroutine: EXPINT [ Ref.5] 


Description: Computes both the matrix exponential 


Fps*ts 
Aps = e 
and the intearal 
Ir Fps*s 
5 е ds 


о 


Reference 5 σ1νθς the method used to perform 
the computations. 


Calling arguements: CALL EXPINT(Fps,Nfps,Aps,Naps,Duml, 
МӘЛІ сс σος Pun?) 
Input arguments: 


Fps Airframe plus actuator plus sensor transfer 
matrix. Generated in ARCRFT 


Nfps Two dimensional row and column vector 
ts Samplina time 
NOD Print parameter: O - Do not print results 


ος ος 
Otherwise print Fps, Ts, е 


Dum2 Vector of working space for computations 


EStput arguments: 


Aps Matrix exponential described above 
Naps Two dimensional row and column vector 
Рим] The integral of the matrix exponential. To be used 


іп computing the discrete control matrix Bps. 


43x10 43x43 43x10 
Bps = Duml x Gps 
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The following subroutines perform the matrix operations 
necessary to compose the A/F-18 system matricies. 


Subroutine: ΠΠ ΟΙ.» N 
Description: Performs the matrix multiplications. 
Calling arguments: CALL MULT(A,NA,B,NB,C, NG) 


input arguments: 
A, B Matricies to be потер ес 


Na, Nb Two dimensional vector giving the number of rows 
and columns in A; “and 72= 


Output arguments: 
С Product of A and B; “С = ЫЕ 


Мс Row and column vector 


Subroutine ADD and SURT are called just as MULT. 
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Eubroutine: NULL [Ref. 5] 


Description:  Generates a null matrix 


Calling arguments: CALL NULL(A,NA) 


input arguments: 
Na Two dimensional vevtor giving the number of rows 
and columns "Of the null matrix, 
Na(l) = Nümber of rows 
Na ο. рее Of columns 


Output arguments: 


A Null matrix having dimension Na(1) x Na(2) 


Subroutine: JUXTR [Ref. 5] 


Description:  Juxtaposes two matricies by row 


Calling arguements: CALL JUXTR(A,Na,B,Nb,C,Nc) 


υπ ὃς arguements: 
А, В The two matricies to Бе juxtaposed by row 


Na, Nb Row and column vectors 


Output arguments: 


А 
С The juxtaposition of A and B С = --- 

В 
Мс Row and column vector 
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Subroutine: JUXT στο πι 


Description:  Juxtaposes two matricies by column 


Calling arguments: CALL JUXTR(A,Na,B,Nb,C,Nc) 


Input arguements: 
A, B The two matricies to be juxtaposed by column 


Na, Nb Row and column vectors 


Output arguments: 


II 
> 
τ 


C The juxtaposition of A and B С 


Νο Row and column vector 


Subroutines: ADD and MULT 
The calling sequences for these subroutines are exactly as 


described for those matrix operations above. 


во 


Eubroutine: OUTPUT 


Bescription: Outputs a matrix with rows and columns 


numbered. 
Шш bing Arguments: CALL OUTPUT(A,Na(1),Na(2),°A ^) 


Pout arguments: 
A Matrix to he output 


Natl), Na(2) Number of rows and number of columns 
in A respectfully 


“A ΄ Matrix name. Must be placed in quotes 
and occupy four spaces. 


Output arguments: 


ШИЕ АҺ matrix 15 written to the FA18 RESULTS file 


lo í 


APPENDIX II 


F/A-18 EXEC PROGRAM 


&TRACE OFF 

CLRSCRN 

GLOBAL TXTLIB VALTLIB VFORTLIB IMSLDP NONIMSL  CMSLIB 
GLOBAL LOADLIB VFLODLIB 

FILEDEF 06 TERMINAL 

FILEDEF Ol DISK Ғ/А-18 DATA Al 

FILEDEF 02 DISK F/A-18 RESULTS (LRECL 132 
FILEDEF 03 DISK OPTMATD DATA Al 

FILEDEF 04 DISK OPTPLOT DA TANAL 

&TYPE LOADING FA18 SIMULATION 

LOAD FA18 ORACLS1 (START 

&EXIT 
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A. 


B, 


ХО 
xW 
ΧΩ 
XW 
XD 
XD 
XD 
XD 
XD 


Е. 


APPENDIX J 


ВЕСНІ CONDITIONS AND STABILITY AND CONTROL DERIVATIVES 


FLIGHT CONDITIONS 


MACH = 0.6 
ALPHA = 2.6184 
PCTHLC = 20.496 
DLHTD = 0.68914 
DLAD = 0.0 
Peer, = 3.0968 
DTEFL = 3.2647 
DRUDL = 0.0 
FLARUD = 0.0 

Бе - 0,2206 
у 0. 1220Е-+Об 


АБТ 


THRUST 
РСТНКС 


DRHD 
DRAD 


DLEFR 
ЕЕЕ ΞΕ 
DRUDR 


DELS 
WAIT 
AX IZ 


T 


B 


H Il 


LONGITUDINAL DERIVATIVES 


ZU 
ZW 


710000.0 
5 
20.496 

0.68914 

ο. ο 
3209068 
3.2647 
ο Ὁ 
ο. 

32550 

ο ο SF ROG 


uH H Hl H 


70 
.196801FE8-01 
.64114Е-02 
-0.45250Е-01 


= 


О 
О 
О 
О 


DATERAL-DIRECTIONAL DERIVATIVES 


т -0,13257Е-01 

= 0.71265E-01 
1). 32650 

КЕ - 0,59729Е-03 
SB = 0.0 
ТН - 0.14257 
S = -0.43034E-03 
LF = -0.412942E-03 
TF = 0.17674Е-03 
= -0.2437 

= 0.0 

= 0.77567 

= 0.67478E-02 
Б << 0.86683Е-03 

= 0.52918E-02 
ШЕ = 0.0 

НТ - -0.21027Е-02 
Ее 0.0 


-0.29125E-01 
= 0.0 

ОШ ООС 
23-1022] 
-0.43808E-02 
0. 12058Е-03 

= 070 

= 0.39478E-02 

= 0.24948 
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ҮТЕСЕ - 646.42 
САМА - 0.0 
Dena 72722770 
XIXZ = -1827.0 
MU = -0.12988E-04 
MW = -0.11331Е-01 
MO = -0.59346 
MWD = -0.34049EF-03 
MDSB = 0.0 
MOTH = 0.0 
MDS - -0.24151Е-02 
МОЕ = -0.5230E-03 
ВО. ЕО 
NV = 0.69481E-02 
NVD = 0.0 
NR = -0.21705 
NP = -0.15254E-01 


NDA = 0.77305E-04 
NDR = -0.46503E-03 
ЕАО. О 

ΠΕΙ 90 .5/361 8-04 
NDTF = -0.32655Е-02 


YARIABLE 


XMACH 
ALT 
VTRUE 
ALFA 
THRUST 
GAMA 
PCTHLC 
PCTHRC 
DLHTD 
DRHTD 
DLAD 
DRAD 
DLEFL 
DLEFR 
ОТЕ 
DTEFR 
DRUDL 
DRUDR 
FLARUD 
DELS8 


WAIT 
XIXX 
ΧΙΥΥ 
Х122 
ХІХ2 


ЕСОН СООО ЕВ 


DEFINITION 


MACH NUMBER 

ALTITUDE, MSL 

TRUE AIRSPEED 

ANGLE OF ATTACK 

TOTAL ENGINE THRUST 

FLIGHT PATH ANGLE 

LEFT THROTTLE POSITION COMMAND (1O04=MAX AB) 
RIGHT THROTTLE POSITION COMMAND(TOOR=MAX АВ) 
LEFT HORIZCNTAL STAB. DEFLECTION (+ TE DOWN) 
RIGHT HORIZONTAL STAB. DEFLECTION (* TE DOWN) 
LEFT AILERON DEFLECTION (+TE DOWN) 

RIGHT AILERON DEFLECTION («ТЕ ООММ) 

LEFT LEADING EDGE FLAP DEFL (+LE DOWN) 

RIGHT LEADING EDGE FLAP DEFL. (+LE DOWN) 

LEFT TRAILING EDGE FLAP DEFL. (+TE DOWN) 
RIGHT TRAILING EDGE FLAP DEFL (+TE DOWN) 

LEFT RUDDER DEFLECTION ( + TE LEFT) 

RIGHT RUDDCR DEFLECTION (+TE LEFT) 

FLARED RUDDER ( PA/TO ONLY ) 

SPEEDBRAKE DEFLECTION 

LONG. CENTER OF GRAVITY LOCATION IN PCT MAC 
AIRCRAFT WEIGHT 

MOMENT Cr INERTIA ABOUT X-BOOY AXIS 

MOMENT OF INERTIA ABOUT Y-BODY AXIS 
MOMENT OF INERTIA ABOUT 2-BODY AXIS 

XZ PLANE PRODUCT Cf INERTIA 
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UNITS 


ND 

FT 
FT/SEC 
DEG 

LB 

RAD 
PERCENT 
PERCENT 
DEG 

DEG 

DEG 

DEG 

DEG 

DEG 

DEG 

DEG 

DEG 

DEG 

DEG 

DEG 
PERCENT 
LB 
SLUG-HIZ 
SLUG IT 
SLUG-FT2 
ΟΠΠ. 


FLIGHT CONDITIONS AND SAMPLING TIME 


AIRCRAFT 


E+05 ALPHA = 0.2618E+01 


Ü 
0 


0Е%00 ALT 
+01 Т5 = 


00 
ОЕ 


OO 


AND FAILURE PARAMETERS 


CONTROL PARAMETERS, 


Sx 


o 


AMPLITUD 


SLARI PY 
RIGHT 


APPENDIX К 


Е ПРЕ ПЕТЬ FILE 
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MATRIX 


GX 


DSD өзе С 
τος 
| 5 
2323 2123 


осоо 


о ΤΟΝΟ 
eS +m 


OOOG 


СЧС 
So oS 
+} i+ 
кїз BR ы 


могло 


OMNGS 
оошо 
a * ? e 
> 


Е 
C S 
+ l| | + 
[13 ΓΤ 3 


— оме о 


СОСО 

oS 
гов 

ообо 
| | 


=“ 


MATRIX 


HX 


ооо 
осо 
+++ 
[τη [ο] [ο] 
σος 
с=с 
e ec 


ооо 


-ro 
SGH 
+ + + 
[τὴ [στ] [1] 
(OCC 
оосо 
=O 


ооо 


oe CN 
= 
Ф + | 
[τ] 153 [τ] 
ΠΟ 
соте 
о т 
4 714% 
ος 
oeg 
Occ 
+ | + 
fr) £523 E 
""ООООО 


oro 


Є © 


NOOO 


ONS 
OOO 


NOM © 


MATRIX 


FYZ 


С) 
ес» 
+ + + + 
[23 23 29 21 
οφ CO OO 
NOOO 
MSc 
аз 9 9 
CIUS eS 


CON e c е 
€ G 
Ф |+ + 
£3 £3 fae 
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